4/QAJT-j 


- /tr - 3£ 

a / 

/d* f - 


FINAL REPORT 
TO 

THE NATIONAL AERONATICS AND SPACE AGENCY 

FOR 

HOLOGRAPHIC GRATINGS FOR SPECTOROGRAPHIC APPLICATIONS: 
STUDY OF ABBERATIONS 

(CONTRACT NO. NGR-2 1-027-010) 


MANJIT S. BHATIA 
PRINCIPLE INVESTIGATOR 


BOWIE STATE COLLEGE 
BOWIE, MARYLAND 

JANUARY 1987 

(NASA-CB- 180267) HCLGGBAPHIC GEA1INGS FCB 
SPECTEGGBAPHIC APPIICATICNS: £1U£I CP 
AEBEBATICHS Final Report (Bo*i€ State 


Cell-, Md.) 


1C 


CSCL 2 OF 


N87-1 5665 


Unclas 

43520 


G3/35 


This is the final report of the project entitled "Holographic 
Gratings for Spectrographic Applications: Study of Abberations" 

conducted at Bowie State College with the help of faculty and 
students. 

This report is divided into the following independent parts, 
each of which can be studied separately: 

1 . Computer Program for Designing Holographic Gratings for 

Seya-Namioka Type Monochromators: Part I: Minimization 

of Abberations Over a Desired Wavelength Range. 

2. Computer Program for Designing Holographic Gratings for 

Seya-Namioka Type Monochromators: Part I: Correcting 

Abberations at Specified Wavelengths. 

3. Design of Astigmatism or Coma Corrected Holographic 
Gratings for Rowland Circle Spectrographic Mounts. 

4. Simultaneous Correction of Astigmatism. Coma, and First 
Spherical Abberation Term of Holographic Concave Gratings 
for Rowland Circle Spectrographs. 

5. Design of Abberation Corrected Holographic Toroid Gratings 
for Seya-Namioka Type Monochromators. 

This report concludes this project. 
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COMPUTER PROGRAM FOR DESIGNING HOLOGRAPHIC GRATINGS 
FOR SEYA-NAMIOKA TYPE MONOCHROMATORS 


INTRODUCTION 

A computer program has been developed based on the theory presented in papers 
by Greiner and Schaeffer (Reference 1) and Noda-Namioka and Seya (Reference 2), 
that determines the optimum holographic grating recording parameters as a 
function of optical instrumental parameters and desired wavelength range for 
Seya-Namioka monochromator mountings. The report is divided into two parts. 

In Part I, the holographic grating recording parameters are determined such 
that astigmatism or coma may be minimized over a selected wavelength range, 
and in Part II, we determine optimum recording parameters such that astigma- 
tism or coma may be corrected to zero at specified wavelengths within the 
desired wavelength range. The program also simultaneously displays the per- 
formance of the holographic grating and equivalent conventional grating as a 
function of wavelength relative to optical aberrations of astigmatism, coma, and 
spherical aberration. 

In this report two methods of designing holographic gratings for operation in 
negative orders only have been developed. The first method is the GENERAL 
METHOD, and with this method the instrumental constants are related only to 
the holographic grating. The second method is called the MODIFIED METHOD, 
and with this method it is possible to design holographic gratings that minimize 
astigmatism or coma over a desired wavelength range for Seya-Namioka mono- 
chromators that are interchangeable with conventional gratings. This means 
that the instrumental constants are the same for both holographic gratings and 
conventional gratings. 

The contents of this report are divided into four sections. Section I presents 
the theory and basic equations for designing holographic concave gratings such 
that minimization of aberrations of astigmatism or coma over a desired wave- 
length range may be achieved. The basic equations are extracted from Refer- 
ence 2 and are presented so that a better understanding of the computer program 
may be appreciated. For more detailed theoretical analysis see Reference 2. 

A description of the computer program including all of the subroutines is pre- 
sented in Section 2. Section 3 presents the input data, definition of terms and 
data card description. Section 4 presents an example for using the program 
and describes the output data in detail. 
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SECTION 1. THEORY AND BASIC EQUATIONS FOR DESIGNING HOLOGRAPHIC 
CONCAVE GRATINGS 


Since a detailed theoretical explanation relative to the design of holographic con- 
cave gratings is given in Reference 2, we only present a brief description for 
ease of understanding the development of the computer program presented in 
this report. 


Figure 1 shows a schematic diagram of the optical system which defines a rec- 
tangular coordinate system. We have used the same terminology as Noda, 
Namioka, and Seya (Reference 2) in order to avoid any ambiguity. Let the ori- 
gin be at vertex O of the concave grating having a radius of curvature R. Let 
the x-axis be the normal to the grating at O and let the xy plane be defined by O 
and two coherent point sources C and D used to record the interference fringes 
on a concave substrate. Points A, P, and B are self-luminous points on the 
entrance slit; a point on the grating and a point at the focus of the diffracted 
image from P of wavelength A in the m-th order, respectively. 


For the ray APB, the light path function F is given by 

F - F ooo + wF ioo + LF on + T w 2 F 2oo 

+ T l2f 020 '•’T^Fjco + twL ! F 1!0 +wLF m 
+ ; w 4 F 4 oo + j w 2 L 2 F, m + i L* F 040 + i wL 2 F M2 

+ 4 F F 022 + T F 03 1 + 7 FF 2 1 1 + * * * 


and 


F ,jk “ M ,jk + * ~ 


/ niX \ 


a) 

( 2 ) 


The M terms refer to the contribution of the conventionally ruled grating and the 
H terms refer to the contribution of the holographic grating. In equation (2), m 
is the order of the grating, X 0 is the laser radiation wavelength, and X is the 
wavelength at which the grating is being employed. The terms M^ k and H jJk 
are defined in terms of the functions f jjk as follows; 

M ijk ss fy k (P.a) + fyi c (p'^). 

^ijk = ^jk^C’^-^jk^D’^ (3) 


for (ijk) = (200), (020), (300), (120), 
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Figure 1. Schematic diagram showing the geometry for recording the holo- 
graphic grating. The instrumental and recording parameters are 
defined as follows : 

Instrument Parameters Recording Parameters 

iy is the object distance from the entrance r c is the distance of recording point source C to 

slit to the center of the grating. the center of the grating. 

iy is the image distance from the exit slit to r D is the distance of the recording point source D 

tite center of the grating. to the center of the grating. 

a is the angle of incident radiation relative 6 is the recording angle for positioning point 

to the normal of the grating. source D relative to the normal of the grating. 

p is the diffracted image angle relative to y is the recording angle for positioning point 

the normal of the giating. source C relative to the normal of the grating. 

A is the entrance slit position and B is the 
exit slit position. 
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( 4 ) 


N, ijk = M 020 -(-'f' n'ijk^.a) + i\ jk (p\0)l, 
H ijk = H 020 -(-I)** l‘ ijk (Pc>7)-l' ljk (P D ,6)l 

fox- (ijk) = (400), (220), (040), where 



r 200^’“) = 

(p cos a - 1) cos a, 





f 020 (P.«> = 

p — cos 0, 





1! 

3 

O 

O 

rO 

p sin a • f 2 oo0°> a )> 





f,2 0 (P> a ) = 

P sin a • f 0 2 0 (P> a )> 





* 400 ^ ,0< ^ " 

P f 200^’ a) l f 200^’ a) ~ 4 P sin2a l> 




f 220 (p ’ a) = 

pf 0 20^> a) sin2 a ~ r 200 

(P.a)], 




f 040 (P.«) = 

Plf 0 20^“)l 2 



( 5 ) 

P = R/r, 

P' = 

; R/r', P c = R/r c , 


R/ r D > 

(6) 

x = r cos 0, y = 

r sin a, x' = r' cos 0, 

y' = 

r' sin 0, 


x c r c 

cos 7, y c 1 

= r c sin 7, x D = r D cos 5, 

1! 

a 

r Q sin 5 . 

( 7 ) 


For equations (2)-(7), R is the radius of curvature of the grating and (r c ,7.0) 
and (r D , 6 , 0) are respectively the cylindrical coordinates of the point sources 
C and D. 

The grating equation is given by 

o (sin oc Q + sin fi 0 ) = mX (8) 

for the principal ray, the ray originating from the center of the entrance slit, and 
diffracted from O at an angle 0. a is the effective grating constant defined by 

1 \) 

a %7rr7„, • < 9 > 

where 6 > 7 . Note that the number of grooves per unit length is N = 1/a. 
Minimization of Aberration Terms 

The theory described up to this point has been presented for the general case. 

We now apply the general theory to a specific application: the Seya-Namioka 
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monochromator. For the Seya-Namioka monochromator the basic instrumental 
parameters are defined as follows (Reference 2); see Figure 2. 


r 0 ,ry = constant, 2k = a Q -0 O = constant 
a 0 = k + 0, (J Q — 0 — k 


( 10 ) 


where r Q is the distance from the center of the entrance slit to O and r^ is the 
distance from O to the center of the exit slit. 2k is the angle AOB and 0 is the 
angle of grating rotation measured from the bisector of the angle 2k and has the 
same sign as the spectral order m. Under conditions (10), the relation between 
0 and X is given by 


X = 



cos k sin 0 


( 11 ) 


We denote F ijk for the ray AOB by F^ . Then, the aberrations in the Seya 
Namioka monochromator can be reduced by minimizing the functions 


F ijk (p.p'.O.k, A iik ) = M iilp + 


ijk - 


ijk 


k ) H * 


" M ijk + A ij|c (sin a 0 + sin 0 O ) 


( 12 ) 


where i + j + k > 2 over a predetermined scanning range, 0 { <0 < 0 2 or 
X, < X < X 2 where 


P = 




(13) 



(14) 


and 0, , and 0 2 are related to X t and X, through equation (11). This is equivalent 
to imposing on F jjk *s the condition 

l ljk = J' 1 F ( 2 ik = minimum (15) 
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Figure 2. Schematic Diagram of the Optical System 
for Seya-Namioka Monochromator 

Determination of Instrumental Parameters p, p', k and A 20 o 

The optimum instrumental parameters are determined in such a way that equa- 
tion (15) is satisfied for the horizontal focusing condition F 200 . This is accom- 
plished by solving the following equations simultaneously, 

^200 o ^*200 

3p 3p 

ai 200 _ ai 200 _ 

dK ’ 3A 200 ’ (16) 
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Once the optimum values of p,p' and x are determined, then each one of the 
integrals I ijk with i + j + k > 2 becomes a function only of A ljk . Then, those 
values of A j j k 's which satisfy equation (15) are calculated from the equation 


where 


^ijk [B ijk (fl 2 )-B ijk (0, )] / [2(0 2 - 0 { ) — (sin 20, “ sin 20 i )1 cos x (17) 

where 

i>2o O (0) = (’P/3)[cos 0 + cos (0 + x) cos (20 + k) 

+ 3 sin (0 + x) sin x] 

+ (p'/3)[cos 0 + cos (0 - x) cos (20 —k) 

— 3 sin (0 — k) sin k] — cos 20 cos k, (18) 

B 0 20^^ = “(P + p ) COS 0 — COS 20 COS K, (19) 


B 3 o O (0) = (p/4 ) 2 (sin (40 + 3k) —4 cos 2(0 + k) sin k — 40 cos x] 

+ (p'/4 ) 2 (sin (40 — 3k) + 4 cos (0 — k) sin k —40 cos k] 

+ (p/6)(3 sin (0 + 2k) — sin (30 + 2x)J 

+ (p'/6)[3 sin (0 — 2k) - sin (30 — 2x)] , (20) 

Bj 20 < 0 ) = TP 2 [sin (20 +x )-20 cosx] 

+ 3 p ' 2 [sin (20 - k) — 20 cos x] 

+ (p/6)(3 sin (0 + 2x) — sin (30 + 2x)] 

( 21 ) 

+ (p‘/ 6)[3 sin (0 - 2 x) - sin (30 - 2 x)] , 

and so on. 
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Determination of the Recording Parameters p c , p D , 7 and 5 and Minimization 
oi' Astigmatism or Coma 

The recording parameters are dependent upon the aberrations to be minimized. 
With a predetermined effective grating constant, the following equations are 
solved simultaneously for minimizing astigmatism and one coma-type aberration, 


sin 6 - sin 7 = \ 0 /a, 

* 200 ^ 0 ^ ~ ^200 - ^ 0 ^ ff ^ 200 ’ 

*020 > 7 ) ~ *020^D - ^0^^020’ 

*300^ 3 C’”^ — ^300^*D = ^0^°^300 ’ (22) 

and for minimizing coma-type aberrations 

sin 5 - sin 7 = X 0 /o, 

*200 ~ * 200 ^* 1 ) = ^ 0 ^°^ 200 ’ 

*300 — *300^D “ ^0^°^300’ 

* 120^0 ~ * 120 ^D ~ ^ 0 ^°^ 120 ' (23) 


Before making any attempt to solve equation (22) or (23) it is necessary to in- 
vestigate the conditions under which equation (22) or (23) can have real solutions 
for p c and p D because H 300 and H 120 are quadratic in p c and p D and their values 
depend oh A 300 and A,, 0 and therefore on values of 0 2 hi I 300 ** 120 * The con- 

dition may be stated in such a way that two quadratic equations of p c or p D re- 
sulting from equations (22) and (23) should not have imaginary roots. To fulfill 
this condition, A 300 and A, 1(| must satisfy the equation 

b 2 — 4 ac > 0 ( 24 ) 


where for A 300 


a = (sin 5 cos 2 7 — sin 7 cos 2 6) cos 2 y, 
b = [ 2 p cos 7 sin 5 + sin (6 — 7) cos 6] cos 7, 
c = pip + cos 6) sin 5 + A 300 cos 2 5 (sin 5 — sin 7), 


(25) 
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and for A 


120 


a = cos 4 y sin 5 — sin y cos 4 6 , 

b = [2 p cos y sin 5 + sin (5 — 7 ) cos 3 6 ] cos 7 , 

c = p (p sin 6 + cos 3 6 ) sin 6 + Aj , 0 cos 4 6 (sin 6 — sin 7 ); (26) 

p = A 200 (sin 6 - sin 7 ) + cos 7 — cos 6 . (27) 

When condition (24) is satisfied, then p c and p D are solutions to equations (22) 
and (23) 

Up to this point, we have been describing the general method of designing aber- 
ration corrected holographic gratins for Seya-Namioka monochromators which 
are not interchangeable with conventional gratings. 


Modified Method 

The modified method of designing holographic gratings for the Seya-Namioka 
monochromator has a practical advantage in that the holographic grating is 
interchangeable with a conventionally ruled grating having the same radius of 
curvature. This is possible provided the same instrumental constants p, p' 
and k used for the conventionally ruled grating are used to design the holographic 
grating. The instrumental constants for the conventional grating are determined 
by solving the equation 

, _ L‘ — , 

' *">00 = J ^200 = minimum 

6 i 

or 

^*200 ^*200 

— — = 0, — r— = 0 

Op dp 

and 

^200 

~k ° (28) 


The design procedure given by equations (15)-(23) must be modified in part in 
order to accommodate condition (28). The modification required is to replace 
equation (16) with equation (28) and equation (17) with (ijk) = ( 200 ) and the rest 
ot the procedure remains the same. 
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SECTION 2 


Section 1 described the theory of designing holographic gratings that minimize 
aberrations over a selected wavelength range. In this section, we use the 
theory and present the computer program which consists of a main program and 
four subroutines. Figure 3 is a diagram of the flow chart for the program. 

PROGRAM DESCRIPTION 

The following is a presentation of the main program and algorithms for each 
subroutine in the computer program to show how the algorithms perform their 
specific tasks. 


MAIN PROGRAM 

The purpose of this program is to direct the flow of calculations depending upon 
which option is selected. If option LMN = 2 is selected, then the optimum value 
of the instrumental parameter angle 2k will be calculated. However, if the 
program user chooses to specify an angle 2k as an input value then option 
LMN = 1 is used. 

The computer program has been developed to accommodate two holographic grat- 
ing design methods; the general method and the modified method. If the general 
method is selected, then the option called IETA = 1 is used in the program. 

When using option IETA = 1, all the instrumental parameters are related only to 
the holographic grating. 

If the modified method is selected, then option IETA = 0 is used in the program. 
When using this option, all instrumental parameters calculated are the same for 
both holographic and conventional gratings. This means that the holographic and 
conventional gratings are interchangeable if they have the same groove frequency 
and radius of curvature. The only difference is that the holographic grating 
minimizes aberrations and the conventional grating does not. 


Once the design method has been selected and the LMN option has been chosen, 
the following steps of calculation are performed. 

STEP A: Read input cards. 

STEP B: Call INSTRM to calculate the optimum angle 2k, and then calcu- 
late Bjj k and A i j k (equation (17) (LMN = 2) ) or calculate the 
functions B^ and A,j k using the input value of 2k (LMN = 1). 
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CALL INSTRM 

CALCULATE INSTRUMENTAL 


PARAMETERS p,p\C 


CALL ASTGMS 
CALCULATE RECORDING 
PARAMETERS 7,fi,P c .P D 


CALL PRFORM 
CALCULATE AND PRINT 
ABERRATIONS (F jjk 's> 


CALL COMA 

CALCULATE RECORDING 
PARAMETERS I^.Pq.Pq 


CALL PRFORM 
CALCULATE AND PRINT 
ALL ABERRATIONS (F jjk 's) 


Figilre 3. Main Program Flow Chart 
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I 

STEP C: Call ASTGMS to calculate the recording parameters p c , p D » 7. j 

and 6 to minimize astigmatism ijk = (020) and one coma term 
ijk = (300). | 

STEP D: Call COMA to calculate the recording parameters p c , P d » 7 > 

and 6 to minimize coma aberrations (ijk) = (300) and (ijk) = (120). 

STEP E: Call PRFORM to display the numerical values of each aberration 
of the holographic grating and compare with that of the equivalent 
conventional grating as a function of wavelength. Note that the 
program automatically calculates the recording parameters for 
both the astigmatism minimization and coma minimization. This 
feature was included in the program so that the program user 
may select which aberration minimization will best produce the 
optimum minimization for all aberrations (astigmatism, coma 
and spherical) throughout the desired wavelength range. 

Subroutine INSTRM 

1 

This subroutine obtains the optimum values of the instrumental parameters, 
namely p, p', 2k and A 200 by simultaneous solution of equation (16). A modified 
form of Newton's method is used to obtain the solution. The quantities B i - k are 
then calculated using the optimum values of p, p', and 2k which in turn give A k 
by equation (17). If a prescribed input value of 2k were used, the selection of 
option LMN = 1 would by-pass the optimization procedures when calculating the 
other instrumental parameters. 

1 

The values of the instrumental parameters as well as A 300 , A 020 , and A 120 are 
printed at the end of this subroutine. 

Subroutine ASTGMS 

The purpose of this subroutine is to obtain the recording parameters such that 
astigmatism (ijk) = (020) and a coma type term (ijk) = (300) are minimized over 
the given wavelength range Xj < X < X 2 . This is achieved by finding a simul- 
taneous solution of equation (22). 

It is first verified that the quantity NX 0 is less than the maximum value of 
I (sin 5 - sin 7)1 , namely 2. If not, an error message is printed stating "error, 

GAMA out of bounds. " 

The values of a, b, c of equation (25) or (26) are then computed and the discrimi- 
nant is checked to ensure b 2 - 4ac > 0. The solutions of the quadratic equation 
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are computed and are checked for po sitiveness; the calculations show that one of 
the two roots [(-b iV b 2 - 4ac )/2a] is always positive and the other always 
negative. 

i . 1 1 - : > • • 

The value of the angle 7 to the desired accuracy is determined by the iterated 
solution ,of the equation (22) by Newton's method. The resulting values of the 
recording parameters 7 and 5 are printed at the end of this subroutine. 


Subroutine COMA 

The purpose of this subroutine is to obtain the recording parameters such that 
the coma type terms (ijk) = (300) and ( 120 ) are minimized over a given wave- 
length range X, < X < X 2 . This is achieved by finding a simultaneous solution 
of equation (23). 

It is first verified that the quantity NX Q is less than the maximum value of 
I (sin £ - sin 7 ) |, namely 2. If not, an error message is printed stating "error, 
GAMA out of bounds. " 

The values of a, b, c of equation (25) or (26) are then computed and the discrimi- 
nant is checked to ensure that b 2 - 4ac > 0. The solutions of the quadratic 
equation are computed and are ver ified for pos itiveness; the calculations show 
that one of the two roots l(-b +Vb 2 - 4ac )/2a] is always positive and the other 
always negative. 

The value of the angle 7 to the desired accuracy is determined by the iterated 
solution of the equation (22) by Newton's method. The resulting values of the 
recording parameters 7 and 5 are printed at the end of this subroutine. 


Subroutine PR FROM 

This subroutine calculates the performance of the holographic grating with the 
recording parameters as computed above and displays the various aberrations 
as a f; motion of the wavelength in the form of a table. The performance of a 
conventional (mechanically ruled) grating with the same groove frequency is also 
shown in the same table for comparison. 

SUCTION 3. DATA SET AND DESCRIPTION OF DATA CARDS 

This section discusses the input data required for a successful operation of the 
program. 
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(a) DATA SET 

The following quantities are required as an input for this program : 

1. \ Q = the wavelength of recording laser light. 

2. X t = the lower limit of wavelength. 

3. X 2 = the upper limit of wavelength. 

4. AX = the wavelength interval selected for displaying the performance 

of the grating as a function of wavelength. 

5. C = the angle 2 k (value 0. 0 is used with option LMN =2). 

6. N = the number of lines per mm (a = 1/N, the grating spacing 

in mm). 

7. R = radius of curvature of the grating. 

8. m = order of the diffraction where m takes on negative values for 

the program, m = -1, -2, etc. 

9. IETA = parameter to specify the grating design method. (IETA - 1) 

specifies the General Method. (IETA = 0) specifies the Modi- 
fied Method. 

10. LMN = parameter to specify whether the angle 2k is to be fed in as 

input (LMN = 1) or is to be computed by the optimization pro- 
cedure (LMN = 2). 

It should be noted "LMN" and "IETA" must be specified in the input for each 
data set for a successful run. 

(b) DESCRIPTION OF DATA CARDS AND SAMPLE DATA CARD 

The data set consists of ten parameters described in (a), The first seven param- 
eters are to be punched in the order shown in (a), starting in column 1 in F10. 0 
format (i. e. , ten columns are reserved for each of the six parameters and each 
of them must include a decimal). 

The parameters m (-1, -2, -3, . . . ), IETA (0,1) and LMN (1,2) take only 
single digit integer values. Two columns are therefore reserved for each of 
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(4) Operating Range = 0 A to 4000 A. 

(5) Order m = -1. 

(6) Determine the optimum angle 2k . 

Solution: (1) indicates that the modified method is used; IETA = 0. (6) indicates 

that the option LMN = 2 is to be used. 

A data card is punched using the data supplied above and punched in the order 
listed below: 

\ 0 = 4579. 3 

X 1 =0.0 

X 2 = 4000. 0 

AX = 500.0 

C = 0. 0 (0. 0 is punched on the input card when the optimum angle 2k 
is to be determined) 

N = 600.0 

R = 1000. 0 

m = -1 

IETA = 0 

, LMN = 2 

Note that the wavelength input values are in angstroms, N in 1/mm, and R in 
millimeters. The angle 2k , where 2k is C on the data card, is in degrees in 
decimal form. 

The output data for the example cited is shown in Figure 5 and described as 
follows : 

Line 1 shows the input parameters; order = -1 is the order of the grating, 

N = groove frequency in grooves per mm, R = radius of curvature 
of the grating in mm, IETA = 0 identifies the design method and 
LMN = 2 shows that the optimum angle 2k has been calculated. 
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Line 2 shows the input values for X 0 = LAMDAO (laser wavelength), 

Xj = LAMDA1 (the lower wavelength limit), and X 2 = LAMDA2 
(upper wavelength limit). 

Line 3 is self-explanatory. 

Line 4 shows the calculated instrumental parameters; C = 2k , p = RHOA, 

P* = RHOB, Tq = SMLR in mm, and rj) = SMPRP in mm. These 
parameters are defined in the text. 

Line 5 shows the type of aberration that has been minimized: Astigma- 
tism = ASTG. 

Line 6 is self-explanatory and refers to line 7. 

Line 7 shows the optimum recording parameters; y - GAMA, 5 = DLTA, 
p c = RHOC and p D = RHOD. These parameters have also been 
defined in the text. 

Line 8 shows the titles for the table showing the performance of the grat- 
ing as a function of wavelength. The table also shows the compari- 
son of the holographic grating with the conventional grating having 
the same instrumental parameters. 

X = LAMDA is the wavelength at which the grating is evaluated. 

F 20 o x 10_4 is the horizontal focus and astigmatism of the first term. 

F 020 is the vertical focus and astigmatism of the second term. 

F 300 is the first coma term. 

F, 20 is the second coma term. 

F 400 is the first spherical aberration term. 

F 220 is the second spherical aberration term. 

F 040 is the third spherical aberration term. 

Line 9 shows the holographic grating performance (Modified Method), the 
value of the wavelength and the values of the aberrations. 
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Line 10 shows the performance of a conventionally ruled grating (MECH. 
RULED). 

Line 11 shows that the following information refers to minimizing coma 
type aberrations for the same instrumental parameters shown in 
Line 4. 

Line 12 is self-explanatory and refers to Line 13. 

Line 13 shows the recording parameters 7 = GAMA, 6 = DLTA, p c = RHOC 
and p D = RHOD required to fabricate a grating to minimize coma. 

Lines 14, 15, and IS are defined in the same way as Lines 8 , 9, and 10 
above. The values in this table refer to a grating designed to 
minimize coma over the desired wavelength. 

REFERENCES 

1. H. Greiner and E. Schaeffer, Optik 16, 350 (1959) 

2. H. Noda, T. Namioka and M. Seya, J. Optical Society of America 64, 

1043 (1974) 

o. T. Namioka, H. Noda and M. Seya, Science of Light 22, 77 (1973) 


19 


CONTENTS 


Page 

INTRODUCTION 1 

SECTION 1. Theory and Basic Equations for Designing 

Holographic Concave Gratings 2 

Minimization of Aberration Terms 5 

Determination of Instrumental Parameters p, p', K and A 7 

Determination of the Recording Parameters p c , p D , y and S to 

Correct for Astigmatism or Coma at a Specified Wavelength 9 

Modified Method 10 

SECTION 2 11 

PROGRAM DESCRIPTION 11 

MAIN PROGRAM 11 

Subroutine INSTRM 13 

Subroutine ASTGMS 13 

Subroutine COMA 14 

Subroutine PR FORM 14 

SECTION 3. DATA SET AND DESCRIPTION OF DATA CARDS 15 

SECTION 4. EXAMPLE AND DESCRIPTION OF DATA 17 


REFERENCES 


20 


COMPUTER PROGRAM FOR DESIGNING 
HOLOGRAPHIC GRATINGS FOR 
SEYA-NAMIOKA TYPE MONOCHROMATORS 

Correcting Aberrations at 
Specified Wavelengths 


Part II: 


COMPUTER PROGRAM FOR DESIGNING 
HOLOGRAPHIC GRATINGS FOR 
SEYA-NAMIOKA TYPE MONOCHROMATORS 

PART II: CORRECTING ABERRATIONS AT 
SPECIFIED WAVELENGTHS 


INTRODUCTION 

A computer program has been developed based on the theory presented in papers 
by Greiner and Schaeffer (1) and Noda, Namioka, and Seya (2), that determines 
the optimum holographic grating recording parameters as a function of optical 
instrumental parameters and desired wavelength range for Seya-Namioka type 
monochromator mountings. In the first report, Part I, the holographic grating 
recording parameters are determined such that astigmatism or coma may be 
minimized over a selected wavelength range. In this report, Part II, we deter- 
mine optimum recording parameters such that astigmatism or coma may be 
corrected at specified wavelengths. By corrected, we mean that the aberrations, 
astigmatism or coma, are zero at specified wavelengths. The program also 
simultaneously displays the performance of the holographic grating and equiva- 
lent conventional grating as a function of wavelength relative to optical aberrations 
of astigmatism, coma, and spherical aberration. In this report, two methods of 
designing holographic gratings for operation in negative orders only have been 
developed. The first method is the general method, and with this method the in- 
strumental constants are related only to the holographic grating. The second 
method is the modified method and with this method it is possible to design holo- 
graphic gratings that correct astigmatism or coma at specified wavelengths within 
the desired wavelength range for Seya-Namioka monochromators that are inter- 
changeable with conventional gratings. This means that the instrumental con- 
stants are the same for both holographic gratings and conventional gratings. 

The contents of this report is divided into four sections. Section I presents the 
theory and basic equations for designing holographic concave gratings such that 
correction for aberrations of astigmatism or coma at a specified wavelength may 
be achieved. The basic equations are extracted from Ref. 2 and are presented so 
that a better understanding of the computer program may be appreciated. How- 
ever, we have modified the equations for minimizing astigmatism or coma such 
that they will now correct astigmatism or coma at specified wavelengths K*. A 
description of the computer program including all of the subroutines is presented 
in Section 2. Section 3 presents the input data, definition of terms, and data card 
description. Section 4 presents an example for using the program and also de- 
scribes the output data in detail. 
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SECTION I. Theory and Basic Equations for Designing 
‘ Holographic Concave Gratings. 

Since a detailed theoretical explanation relative to the design of holographic con- 
cave gratings is given in Ref. 2 we only present a brief description for ease of 
understanding the development of the computer program presented in this report. 

Figure 1 shows a schematic diagram of the optical system which defines a rec- 
tangular coordinate system. We have used the same terminology as Noda, Namioka, 
and Seya (Ref. 2) in order to avoid any ambiguity. Let the origin be at vertex 0 of 
the concave grating having a radius of curvature R. Let the x-axis be the normal 
to the grating at 0 and let the xy plane be defined by 0 and two coherent point 
sources C and D used to record the interference fringes on a concave substrate. 
Points A, P, and B are self luminous points on the entrance slit, a point on the 
grating and a point at the focus on the diffracted image from P of wavelength 
in the m th order, respectively. 

For the ray APB, the light path function F is given by 

F “ F 000 + wF 100 + tfou + ^ w2F 200 + ~2 ^ F 020 + *2 w2F 300 


4 wliF .» * “ 1F m + j* <F «o *7 * j 1 2f M 0 4 l4F »<° 

4' 1,,F *« 4 llF °« * T l3F »3t 4" ,lF ,u ♦ U) 

and 

F, jk =M iik .(»\A 0 )H lik . (2) 

The M terms refer to the contribution of the conventionally rulled grating and 
the H terms refer to the contribution of the holographic grating. In Eq. 2, m is 
the order of the grating, k 0 is the laser radiation wavelength and k is the wave- 
length at which the grating is being employed. The terms Mijk and Hijk are de- 
fined in terms of the functions fijk as follows: 

M ijW = 

H .jk = f ijk (P c >7)- f ijk (pD,S) (3) 


2 


for 

(ijk) = (200), (020), (300), (120), 

“ iik =M„J0 - t f 1Jk <^./3)] , 

H, jk --H 0!0 -(-l)«Mf lik (pc.r)-f. jk (^D,S)] (4) 

for 

(ijk) = (400), (220), (040), 

where 

f 200 (p,a) = (pcosa-1) cosa, 

^ 020 ( ^ ’ a ) — P “ cosa, 
f 30o(P' a '> = P Sina ‘ f 2 00 (^- a )- 

^120 ’ a ) ~ P S *- na * ^020 (P * a ) * 

f 400^’ a ) = P f 20<)(P’ a ) tf 200 (P' a ) ~ 4 P sir > 2 a] • 

^220 (P ’ a ^ ~ P ^020 (P ' a ) ^ sin 2 a - ^200 • a )^ > 

f 04o(^' a ) = pt f 020<^ ,a )^ 2 < 5 ) 

p = R/r , p' = R/r' . p c = R/r c , p D = R/r D , (6) 

x = r cos a, y': r sina, x' = r'cos/ 3 , y' = r' s in / 3 , 

x c =r c cos7, y c = r c s i n 7 , x D = r Q cos 8 , y D =r D sinS. ( 7 ) 

For equations (2)-(7), R is the radius of curvature of the grating and (r c , y, 0) 
and (r D , 8, 0) are respectively the cylindrical coordinates of the point sources 
C and D. 

The grating equation is given by 

a( sin a Q + sin fi Q ) = m\ (8) 

for the principal ray, the ray originating from the center of the entrance slit, 
and diffracted from 0 at an angle fi. a is the effective grating constant defined by 
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Figure 1. Schematic diagram showing the geometry for recording holographic 
gratings. The instrumental and recording parameters are defined as follows: 

(Instrumental Parameters): 

r 0 is the objective distance from the entrance slit to the center of the grating, 
r' is the image distance from the exit slit to the center of the grating, 
a is the angle of incident radiation relative to the normal of the grating, 
p is the diffracted image angle relative to the normal of the grating. 

A is the entrance slit position and B is the exit slit position. 

(Recording Parameters) : 

r c is the distance of recording point source C to the center of the grating. 
r d is the distance of the recording point source D to the center of the grating. 

5 is the recording angle for positioning point source D relative 
to the normal of the grating. 

> is the recording angle for positioning point source C relative 
to the normal of the grating. 
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a 


sin S - 


s in y 


where 


8 > y 


Note that the number of grooves per unit length is N = l/a. 


Minimization of Aberration Terms 

The theory described up to this point has been presented for the general case. 
We now apply the general theory to a specific application; the Seya-Namioka 
monochromator. For the Seya-Namioka monochromat the basic instrumental 
parameters are defined as follows; (reference 2) see Figure 2. 

r 0> = constant , 2K = a Q - /j Q = constant 

a 0 = K + (9, /S o = 0-K (10) 

where r 0 is the distance from the center of the entrance slit to 0 and rj, is the 
distance from 0 to the center of the exit slit. 2K is the angle AOB and Q is the 
angle of grating rotation measured from the bisector of the angle 2K and has the 
same sign as the spectral order m. Under conditions (10), the relation between 
d and \ is given by 


\ = — cos K sin d (11) 

We denote F ijk for the ray AOB by F ijk . Then, the aberrations in the Seya- 
Namioka monochromator can be reduced by minimizing the functions 


F. 


mX. 


i j k ( * P ' K, A i . k ) = + Hj. k 

% 


= M ijk + A ijk( sina o + sin ^o) ( 12 ) 


where 


i + J + K ^ 2 
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z 



Figure 2. Schematic diagram of the optical system for 
Seya-Namioka monochromator 
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over a predetermined scanning range, d x = 6 i d 2 or = X = X 2 where 



( 13 ) 


H 


ijk 


‘i* sinS-siny 


■ft)"- 


(14) 


and 0, , and 0, are_related to X, and \ 2 through equation (11). This is equiva- 
lent to imposing on F Jjk 's the condition 


ijk 


r 


= minimum 


(15) 


Determination of Instrumental Parameters p, p‘ , K and A 200 


The optimum instrumental parameters are determined in such a way that Eq. (15) 
is satisfied for the horizontal focusing condition F 200 . This is accomplished by 
solving the following equations simultaneously, 


31 go 0 21 mo. 

.if, ' Tip, 


200 _ 0 200 _ Q 

BK = ’ 3A 200 - 


(16) 


Once the optimum values of p, ~p ' and K are determined, then each one of the 
integrals I . with i + J + K = 2 becomes a function only of A t jk . Then, those 
values of A** k 's which satisfy Eq. (15) are calculated from the equation 


where 




[2(d 2 - 0j) - (sin 2^ 2 - sin20j )] cosK, 


(17) 
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where 


B 2OO (0) = (P /3 > [cos0 + c °s(0 + K) cos(20 + K 
+ 3 sin(<9 + K) sinK] 

+ (/?' /3 ) [cosd + cos (0 - K) cos(20 - K) 

- 3 sin(# - K) sinK] - cos20 cosK, (18) 

B 020 ($) = 2 (p + ~p‘ ) cosd - cos2<9 cosK, (19) 

B 300 (P) = (^ 4 ) 2 tsin(40 + 3K) 

- 4 cos2(0 + K) sinK - 4 9 cosK] 

+ (P7 4) 2 [ s i n ( 4(9 - 3K) 

+ 4 cos (0 - K) sinK - 4(9 cosK] 

+ (p/6 ) [3sin(0 + 2K) - sin(3tf + 2K)] 

+ (p'/6) [3 sin(0 - 2K) - sin(3<9 - 2K)] , (20) 

B no( (? ) = |p 2 [sin(2<9 + K) - 2(9 cosK] 

+ -|p' 2 [sin(20 - K) - 29 cosK] 

+ (p/6) [3 sin(0 + 2K) - sin(3P + 2K)] 

+ (p'/6) [3 sin(0 - 2K) -sin(30 - 2K>] , (21) 

and so on. 
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Determination of the Recording Parameters p c » p D , y and s to Correct 
for Astigmatism or Coma at a Specified Wavelength 

The recording parameters are dependent upon the aberrations to be corrected 
and the selected correcting wavelength \* . With a predetermined effective grat- 
ing constant, the following equations are solved simultaneously for correcting 
astigmatism and one coma type aberration, 


sin 8 - sin 7 = NX. 0 

*20o(Pc ' y ) “ *200 (p D ,S) = N\ 0 A 

200 

k Q 

* 30 o(Pc • T) “ * 200 • 8 ) ~Z [p s ina* (p cos 2 a* - cos a* ) 

m K 

+ p' sin/ 3 *( p' cos 2 / 3 * - cos/ 3 * )] 

X 0 

*02o(P c - v) - *020 (Po- 2 ) = - — [p + p' - (cosa* + cos/ 3 *)] 

m\* 


\ 


> 

/ 


( 22 ) 


and for correcting coma-type aberrations 

s in S - sin 7 = NX 0 ^ 

*200 (P c - y ) “ *200 (Pd > 2 ) = N\q ^200 

k 0 

*300 ( Pc * y) - *300 ( Pq ’ 8 ) rr Ip sina*(pcos 2 a* - cosa* ) 

m A.* 

> (23) 

+ p' sin/ 3 *(p' cos 2 / 3 * - cos/ 3 *)] | 

\ 0 

*120 (Pc - y) “ * 12 o(Pd’ 8 ) 71 lpsina*(p - cos a*) 

m K 

+ p' s in/ 3 *(p- - cos/ 3 * )] / 

Before making any attempt to solve Equations (22) or (23) it is necessary to in- 
vestigate the conditions under which Equations (22) or (23) can have real solutions 
forp c and p D because H 300 and H 120 are quadratic in p c andp D mid their val- 
ues depend on A 300 and A 120 and therefore on values of t> 2 in I 300 and I l20 . 

The condition may be stated in such a way that two quadratic equations of p c or 
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P D resulting from equations (22) and (23) should not have imaginary roots. To 
fulfill this condition, A 300 must satisfy the equation 


b 2 - 4ac l 0 

where for A 

300 


(24) 


a = ( s in S cos 2 ')' - sinycos a $) cos 2 y , 


b = [2pcosysinS + sin(S-y)cosS] cos y , 
c = p(p + cos 5 ) s in S + A 200 cos 2 S(sin8 - s iny) , (25) 


p = A 200 (sinS -siny) + cosy-cosS. (26) 

When condition (24) is satisfied, then p c and p Q are solutions to Equations 22 and 
23. Up to this point, we have been describing the general method of designing 
aberration corrected holographic gratings for Seya-Namioka monochromators 
which are not interchangeable with conventional gratings. The next paragraph de- 
scribes the method where the conventional and holographic gratings are inter- 
changeable. 


Modified Method 


The modified method of designing holographic gratings for the Seya-Namioka 
monochromator has a practical advantage in that the holographic grating is inter- 
changeable with a conventionally ruled grating having the same radius of curvature. 


This is possible provided the same instrumental constants p, p> and K used for 
the conventionally ruled grating are used to design the holographic grating. The 
instrumental constants for the conventional grating are determined by solving the 
equation 


or 


200 



dd = minimum. 


and 


^ ^ 200 A ^ ^ 200 
*P = ' Tip 1 



( 27 ) 
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The design procedure given by Equations (15)-(23) must be modified in part in 
order to accommodate condition (27). The modification required is to replace 
Eq. (16) with Equation (27) and Equation (17) with (ijk) = (200) and the rest of 
the procedure remains the same. 


SECTION 2 

Section 1 described the theory for designing holographic gratings that correct 
astigmatism or coma at specified wavelengths within a desired wavelength range. 
In this section, we use the theory and present the computer program which con- 
sists of a main program and four subroutines. Figure (3) is a diagram of the 
flow chart for the program. 


PROGRAM DESCRIPTION 

The following is a presentation of the main program and algorithms for each sub- 
routine in the computer program to show how the algorithms perform their spe- 
cific tasks. 


MAIN PROGRAM 

The purpose of this program is to direct the flow of calculations depending upon 
which option is selected. If option LMN = 2 is selected, then the optimum value 
of the instrumental parameter angle 2K will be calculated. However, if the pro- 
grammer chooses to specify an angle 2K as an input value, then option LMN = 1 
is used. 

The computer program has been developed to accommodate two holographic 
grating design methods; the general method and the modified method. If the 
general method is selected, then the option called IETA = 1 is used in the pro- 
gram. When using option IETA = 1, all the instrumental parameters are related 
only to the holographic grating. 

If the modified method is selected, then option IETA = 0 is used in the program. 
When using this option, all calculated instrumental parameters are the same for 
both holographic and conventional gratings. This means that the holographic and 
conventional gratings are interchangeable if they have the same groove frequency 
and radius of curvature. The only difference is that the holographic grating cor- 
rects aberrations to zero and the conventional grating does not. 
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INPUT: X Q , X,, X 2 , K, AX. C, X*. N, R 
LMN, m, I ETA 



Figure 3. Main Program Flow Chart 


12 








Once the design method has been selected and the LMN option has been chosen, 

the following steps of calculation are performed. 

STEP A: Read input cards. 

STEP B: Call INSTRM to calculate the optimum angle 2K, and then calculate Bijk 
and Aijk Eqs. (17) (LMN 2) or calculate the functions Bijk and Aijk using 
the input value of 2K (LMN = 1). 

STEP C: Call ASTGMS to calculate the recording parameters p c , p D ,y, and 2 
to correct astigmatism ijk = (020) and one coma term ijk = (300) at \* . 

STEP D: Cali COMA to calculate the recording parameters P c , p D >y, 8 to cor- 
rect coma aberrations (ijk) = (300) and (ijk) = (120) at K*. 

STEP E: Call PR FORM to display the numerical values of each aberration for 
the holographic grating and compare with that of the equivalent con- 
ventional grating as a function of wavelength. Note that the program 
automatically calculates the recording parameters for both the astig- 
matism correction and coma correction. This feature was included in 
the program so that the program user may select which aberration 
correction will best produce the optimum performance for all aberra- 
tions (astigmatism, coma and spherical) throughout the desired wave- 
length range. 


Subroutine INSTRM 

This subroutine obtains the optimum values of the instrumental parameters 
namely p, p', 2K and A200 by simultaneous solution of Eqs. (16). A modified 
form of Newton's method is used to obtain the solution. The quantities Bijk's are 
then calculated using the optimum values of p , p ‘ and 2K wliich in turn give A ; t k 
by Eq. (17). If a prescribed input value of 2K were used, the selection of option 
LMN = 1 would by-pass the optimization procedure when calculating the other 
instrumental parameters. 

The values of the instrumental parameters as well as A 300 , A oao , and A 1J0 are 
printed at the end of this subroutine. 


Subroutine ASTGMS 

The purpose of this subroutine is to obtain the recording parameters such that 
astigmatism (ijk) = (020) and a coma type term (ijk) = (300) are corrected for a 
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specified wavelength A.* . This is achieved by finding a simultaneous solution of 
Eq. (22). 

It is first verified that the quantity NA . 0 is less than the maximum value of 
| (sin S - siny) | ; namely 2. If not, an error message is printed stating "error, 
gamma out of bounds." 

The values of a, b, c of Eqs. (25) are then computed and the discriminant is checked 
to ensure b 2 - 4ac = 0 . The solutions of the quadratic equation are computed and 
c hecked fo r positiveness; the calculation shows that one of the two roots [-b + 

Vb 2 - 4ac/2al is always positive and the other always negative. 

The value of the angle y to the desired accuracy is determined by the iterated 
solution of Eq. (22) by Newton's method. The resulting values of the recording 
parameters y and 8 are printed at the end of this subroutine. 


Subroutine COMA 

The purpose of this subroutine is to obtain the recording parameters such that 
the coma type terms (ijk) = (300) and (120) are corrected at a specified wave- 
length K*. 

It is first verified that the quantity N \ 0 is less than the maximum value of 
[ (sin 8 - sin 7 ) |; namely 2. If not, an error message is printed stating "error, 
gamma out of bounds." 

The value of a, b, c of Eq. (25) are then computed and the discriminant is checked 
to ensure that b 3 - 4ac = 0. The solutions of the quadratic equation are computed 
and ver ified for positiveness; the calculation shows that one of the two roots 
[-b + Vb 2 - 4ac/2a] is always positive and the other always negative. 

The value of the angle y to the desired accuracy is determined by the iterated 
solution of Eq. (22) by Newton's method. The resulting values of the recording 
parameters y and 6 are printed at the end of this subroutine. 


Subroutine PRFORM 

This subroutine calculates the performance of the holographic grating with the 
recording parameters as computed above and displays the various aberrations 
as a function of the wavelength in the form of a table. The performance of an 
equivalent conventional (mechanically ruled) grating with the same groove fre- 
quency is also shown in the same table for comparison. 
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SECTION 3. DATA SET AND DESCRIPTION OF DATA CARDS 


This section discusses the input data required for successful operation of the 
program. 


(a) DATA SET 


The following quantities are required as input data for this program: 


(1) 

^0 

(2) 


(3) 

s 

(4) 

A \ 

(5) 

\* 

(6) 

C 

(7) 

N 

(8) 

R 

(8) 

LMN 

(10) 

m 

(11) IETA 


wavelength of laser recording radiation, 
lower wavelength limit, 
upper wavelength limit. 

wavelength interval at which aberrations will be displayed, 
wavelength at which aberrations will be corrected, 
angle 2K (value 0.0 is used with option LMN = 2). 
number of lines per mm. (a= 1/N; grating spacing), 
radius of curvature of the grating in mm. 

Parameter to specify whether the angle 2K is to be fed in as an 
input (LMN = 1) or is to be computed by the optimization procedure 
(LMN = 2). 

order of diffraction where m takes on negative values. 

Parameter to specify the grating design method: 

(IETA - 1) specifies the general method 
(IETA = 0) specifies the modified method. 


Note that "LMN" and "IETA" must be specified in the input for each data set for 
a successful run. 


(b) DESCRIPTION OF DATA CARDS AND SAMPLE DATA CARDS 

The data set consists of eleven parameters described in (a). The first eight 
parameters are punched on the first card in the order stated in (a) starting in 
column 1 in F10.0 format. For example, ten columns are reserved for each of 
the eight parameters and each must include a decimal. The parameters LMN 

(1,2), ra(-l, -2, ) and IETA (0, 1) are punched on the second card using 

only single digit integer values and two columns are reserved for each of these 
parameters starting in column 1. Figure 4 illustrates typical data cards with 
the input values shown at the top. 
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Note that as many data sets as desired may be put into a single program run 
provided that the values of all eleven parameters are punched for each data set 
which consists of two cards. 


SECTION 4 

This section presents a typical example for using the program and also describes 
the output data in detail. 

Example: Determine the optimum instrumental and recording parameters for 
the design of a holographic grating that corrects astigmatism or coma to zero 
at 2500 A having the following specifications: 

(1) Groove frequency = 600 1/mm 

(2) Radius of curvature - 1000 mm 

(3) Operating range = 0 A to 4000 A 

(4) Order m - -1 

(5) Determine the optimum angle 2K 

(6) This holographic grating is not interchangeable with a conventional grating. 

Solution: (5) indicates that the option LMN = 2 will be used. (6) specifies that 
the General Method of design will be used (IETA = 1). A data card is punched 
using the data supplied above and punched in the order listed below: 

k Q = 4579.3 

\j =0.0 

a- 2 = 4000.0 
A \ = 500.0 
\* = 2500.0 

C = 0.0 (0.0 is punched on the input card when the optimum angle 2K is 
to be determined). 

N = 600.0 
R = 1000.0 
LMN = 2 
m = -1 
IETA = 1 

Note that the wavelength input values are in angstroms, N in lines per mm and 
R in mm. The angle 2K, where 2K = C on the data card, is in degrees in decimal 
form. 

The output data for the example cited is shown in Fig. 5 and described as follows: 


17 


Line 1; shows the input parameters; order = -1 is the order of the grating, 

N = groove frequency, R = radius of curvature of the grating, IETA = 1 identi- 
fies the design method, and LMN = 2 shows that the optimum angle 2K has been 
calculated. 

Line 2; the input values for k Q = LAMDAO (laser wavelength), A- t = LAMDA1 
(lower wavelength limit), k 2 ~ LAMDA2 (upper wavelength limit) and k* = LAM- 
DAC (correcting wavelength). 

Line 3; is self explanatory. 

Line 4; shows the calculated instrumental parameters C = 2K, p = RHOA, p‘ = 
RHOB, r^ = SMLR in mm and r Q ‘ = SMPRP in mm. These parameters are de- 
fined in the text. 

Line 5; shows the type of aberration that has been corrected; Astigmatism = ASTG. 
Line 6; is self explanatory and refers to line 7. 

Line 7; shows the optimum recording parameters y = GAMA, & = DLTA, p c = 

RHOC and p Q = RHOD. These parameters have also been defined in the text. 

Line 8; shows the wavelength (k = LAMDA) at which the grating is evaluated and 
the aberration terms described as follows: 

F 200 x 10 -4 is the horizontal focus and astigmatism of the first term. 

F 020 is the vertical focus and astigmatism of the second term. 

F 300 is the first coma term. 

Fi 20 is the second coma term. 

F 400 is the first spherical aberration term. 

Fj 20 is the second spherical aberration term. 

F 040 is the third spherical aberration term. 

Line 9; shows the holographic grating performance (General Method), the evalu- 
ating wavelength and the values of the aberrations. 

Line 10; shows the performance of the equivalent conventional grating (MECH. 
RULED). 
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Line 11; shows that the following information refers to correcting coma type 
aberrations for the same instrumental parameters shown in line 4. 

Line 12; is self explanatory and refers to line 13. 

Line 13; shows the recording parameters 7 = GAMA, S = DLTA, p c = RHOC, 
P D = RHOD required to fabricate a grating to correct for coma. 

Lines 14, 15 and 16 are defined in the same way as lines 8, 9, and 10 above. 
The values in this table refer to a grating designed to correct for coma at a 
specified wavelength. In this example the correcting wavelength is 2500A. 
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DESIGN OF ASTIGMATISM OR COMA CORRECTED HOLOGRAPHIC 
GRATINGS FOR ROWLAND CIRCLE SPECTROGRAPH 1C MOUNTS 


INTRODUCTION 

For many spectroscopic instruments such as spectrographs, monochromators, spectropho- 
tometers, and similar instruments employing concave optical reflecting diffraction gratings, 
a major inherent optical disadvantage is that the diffracted image size (exit slit) is distorted 
and many times larger than the objective size (entrance slit). This phenomenon is ba- 
sically due to astigmatism caused by, in large part, the conventional concave diffraction 
grating used in the instrument. In order to detect all the energy in such large diffrac- 
ted images, very large area detectors must be used. These detectors, if available, would 
be very large in size and would be unacceptable for use especially in flight spectroscopic 
instruments. At the present time, small area detectors are used with much of the energy 
being lost. A solution to this problem is to use the aberration corrected holographic grat- 
ings that correct for astigmatism. 

The purpose of this report is to present the design of holograplric gratings that will in 
many cases solve the problem of astigmatism and in other cases the problem of coma 
which causes a loss in resolving power in spectroscopic instruments. 

The first section of this report presents the theory and basic equations for designing holo- 
graphic diffraction gratings. The remaining sections are devoted to the development of 
a computer program that calculates the necessary geometrical data which is one of many 
requirements that makes the fabrication of holographic gratings a reality. The compu- 
ter program has been developed such that optimum recording parameters as a function 
of instrument parameters are determined that make astigmatism or coma zero at speci- 
fied wavelengths for spectrographs based on the Rowland circle mount. When optimum 
recording parameters have been determined, the program will display, for the purpose 
of comparison, numerical tables of the performance for the holographic grating and 
equivalent conventional grating as a function of wavelength. This numerical comparison 
makes it possible to analyze the data and decide whether or not the holographic grating 
performance is better than the equivalent conventional grating relative to optical aberrations. 


THEORY AND BASIC EQUATIONS FOR DESIGNING HOLOGRAPHIC SPHERICAL 
CONCAVE GRATINGS 


Figure 1 shows the optical system for using and recording holographic gratings. The rec- 
tangular coordinate system is defined as follows: I is the origin at the vertex of the grat- 
ing blank and the x-axis is normal to 1. The x-y plane is defined by I and the two recording 
point sources C and D. Points A, M, and B are respectively, a self luminous point on 
the entrance slit, a point on the grating surface, and a point on the focused diffracted 
image from M at a wavelength X in the n 1 * 1 order. As a result of this geometry the 
aberrant optical path function for a conventional grating may be written as 


1 



z 



Figure 1. Schematic diagram of the optical system. A is a self luminous point on the 
entrance slit, M is a point on the n 1 * 1 groove and B is the spectral image point. C and 
D are the recording point light source positions. The performance and recording param- 
eters are defined as follows: 

Performance Parameters : 

La is the distance from A (entrance slit) to the center of I of the grating. 

Lg is the distance from the diffracted image B to the center 1 of the grating. 
a is the angle of incidence. 

P is the diffraction angle, 
z is' the object height. 
z is the image height. 

Recording Parameters: 

Lc is the distance from the point source C to the center I of the grating. 

Lj) is the distance from the point source D to the center I of the grating. 

7 is the position angle for point source C relative to the normal to the grating. 

6 is the position angle for point source D relative to the normal to the grating. 
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AM = AM + MB - (1A + IB) - KNX 


( 1 ) 


where K is the order of diffraction, X is the diffracted image wavelength and N is the grating 
groove frequency. 

When a grating blank coated with the photoresist is illuminated with coherent spherical waves 
originating from points C and D, grooves are recorded, after developing, in the photo- 
resist in response to the interference fringes produced at the site of the grating blank. 
The interference fringes are formed according to the equation 

N = MC - MD - (1C - IP) (2) 


where X 0 is the recording laser wavelength. 

Substituting eq. (2) into eq. (1) gives the aberrant optical path function 

AM = MA + MB-(IA + IB)-^ j(MC - MD) - (1C - ID) 
for the holographic grating. 


(3) 


The quantities on the right side of AM in eq. (3) are derived based on the coordinate system 
described in Figure 1 and the equations of a sphere. Since the results of these deriva- 
tions are lengthy, we suggest that the reader see reference 1 for more detailed analysis. 
Therefore, it sufficies to say that after a series expansion is performed, the following 
aberrant optical path function AM is obtained. 


AM 


= -Ylsina + sin/3-^ {sin 7 -sin 6 } l+^lAj] 


Y 2 


(4) 


+Y [A 2 ] +Y ICi J +^IC 2 ] +Y4 ISj 1 
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Note that in the expressions of equation (5), L^, Lg, a and /3 are the terms for the conven- 
tionally ruled grating, and Lq, Ljj, 7, and 8 are the terms related to the holographic contri- 
bution. R is the radius of curvature of the grating, K is the order of the diffracted 
wavelength, and X Q is the recording wavelength. The aberration terms are defined as 
follows; Aj is the horizontal diffracted image focus, A2 is the astigmatism or vertical 
focusing term, Cj and C2 are coma terms and Sj , S2, and S3 are spherical aberration terms, 

The conventional grating equation is given by 


a (sin a + sin 0) = KX 


( 6 ) 


and for the holographic grating o is the effective grating constant and is defined by 
equation 


sin 7 - sin 8 


where 7 > 6 


(7) 


The preceding information describes the basic general theory for holographic gratings. What 
follows is an application of the general theory for a specific spectroscopic instrument; cor- 
recting astigmatism or coma at a specified wavelength for a Rowland circle spectrographic 
mount. 


HOLOGRAPHIC GRATING DESIGN FOR THE ROWLAND CIRCLE SPECTROGRAPH 


In order to design a holographic concave grating that has the Rowland circle as its fo- 
cal curve; it is necessary to make the horizontal focusing term Aj and the coma term 
Cj zero at all wavelengths in equation (5). This condition is satisfied by 

L^ = R cos a, Lg = R cos (3, Lq = R cos 7 and Lq = R cos 8 (8) 


The aberrant optical patli function for condition (8) will now have the form 
AM =-Y[sin a + sin | sin 7 - sin 5 J J [A^'l 

YZ‘ fp Y* ,, Y* Z- K ,, 3 

+ 2Rl [ - 1 + 8R3 I S 1 J + W I S 2 ) + srJ l s I J J 


(9) 


and the expressions for astigmatism an< * coma C2* in equation (5) reduce to 


, _ sin^ or* 
- cos a* 


sin^jg* KX* 
cos/J* “ X 0 


sin- 7 
cos 7 



( 10 ) 
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and 


C 2 ' = 


sin-* g* 
cos2 a* 


+ sin^ 0 * 
cos2 0* 


KX* sin^ 7 sin^ 5 

Xo cos2 y cos2 8 


(ID 


where X* is the wavelength at which the aberrations are to be corrected to zero and 0* is 
calculated from 


a (sin a* + sin 0*) - KX* (12) 

The recording parameters y, 8, L c , and L D are determined for a chosen set of values o, X*, 
and a* to eliminate astigmatism A 2 or coma c{ at a specified wavelength X*. If astigma- 
tism A 2 ' is desired to be eliminated at X*, then the following set of equations must be 
solved; 


a (sin a* + sin 0*) = KX*, (12) 


sin 7 - sin 6 v ’ 

and equation (10), A 2 ' = 0. For this case, the recording angles y and 8 are chosen such 
that astigmatism A 2 ' is made zero at X*. 

If one chooses to eliminate coma at a wavelength X*, then the following set of equa- 
tions must be solved; (12), (7) and (11), C 2 ! = 0. For this case, the recording angles 
7 and 5 are chosen such that C 2 ' is made zero at X*. 

When the recording parameters y and 8 have been determined for a given set of a, a*, 
and X*, then the values of L^, Lg, Lc and Lq are calculated from equation (8). The 
remaining aberrations Si',S 2 ', and S3' in equation (5) are calculated after the aberra- 
tions A 2 ' and C 2 ' have been made zero at X* for the values determined for 7 and 6. 
The aberration calculations determined by the computer for astigmatism, coma, and 
spherical are independent of Y, Z, R, and the constants in the denominators of the ex- 
pressions for the coefficients. In other words the values of the aberrations are pure 
numbers. 


OVERALL LOGIC 


In the first section, we presented the general theory of holographic concave gratings and ap- 
plied the theory to designing concave holographic gratings, corrected for astigmatism or 
coma, for a spectrograph based on the Rowland circle mount. 

In this section, we present the logical steps of the computer program that leads to the 
design of holographic gratings. 
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Step 1 requires that one selects the groove frequency, angle of incidence, correcting 
wavelength X*, and the radius of curvature R for the holographic grating. 

Step 2 ; Based on the data of step 1 , determine the angle of diffraction 0 * at X*. 

Step 3 ; The expression for astigmatism A2' equation ( 10 ) is then scanned as a function 
of the recording angle 7 to determine the optimum angles 7 and 6 that will make A2' = 0 
at X*. Successive values of 7 are obtained by an equation of the form 7 = 7j + A7 where 
A7 is an incremental angle and chosen to equal 0 . 25 °. Therefore the accuracy of 7 is de- 
pendent on the value of A7; the smaller the incremental angle Ay the better the accuracy 
of 7. 

Step 4 ; When the optimum recording angles 7 and 6 have been determined for correcting 
astigmatism A2' at X*, then the performance of the grating relative to all aberrations as a 
function of wavelength over a chosen wavelength range is calculated for the holographic 
grating and displayed in the form of a table. For comparison, the performance of an 
equivalent conventional grating is also calculated and displayed for the same wavelength 
range. 

The same logical process is employed for correcting coma C2' (equation ( 1 1 )), however, 
the optimum recording parameters are related to correcting coma. 

A unique feature of the program is that it automatically calculates the optimum record- 
ing parameters for both astigmatism and coma corrections and then displays complete 
data for each aberration correction. This feature was included in the program so that the 
program user may select which aberration correction will best provide the optimum per- 
formance for all aberrations (astigmatism, coma, and spherical) throughout the desired 
wavelength range. 


DESCRIPTION OF COMPUTER PROGRAM 


The previous section gave the logical sequence for setting up the computer program. This 
section describes the computer program which consists of a main progrant and two sub- 
routines. 

The main program directs the How of calculations and calls upon the subroutines to per- 
form their various functions. The flow diagram (Figure 2 ) shows that the following events 
will be performed. 

1 . Input data cards are read. 

2 . Call CORECT. CORECT is a subroutine that determines the optimum recording 
parameter values for 7 and S such that astigmatism is corrected to zero at X*. This 
is achieved by simultaneous solution of equations ( 12 ), ( 7 ), and ( 10 ). Then this sub- 
routine also automatically determines the optimum recording parameter values for 7 
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and 6 to correct coma to zero at X*. This task is accomplished by simultaneous so- 
lution of equations (12), (7), and (11). 

3. Call PRFORM. PRFORM is a subroutine that calculates the numerical values lor all 
aberrations as a function of wavelength for both holographic and conventional gratings. 

4. Output. This function displays the numerical values of the input data, recording 
parameters, and all aberrations for the holographic and conventional gratings as a func- 
tion of wavelength in the form of a table. 

MAIN PROGRAM FLOW DIAGRAM 


START 


I 

INPUT 

R, K, N, X 0 ,X\ X u ,Xu 
a L , a U' Aa 


♦ 





Figure 2. Main Program Flow Diagram. 
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INPUT DATA, DEFINITION OF TERMS AND DATA CARDS 


This section describes the input data and defines the parameters that will be punched 
on the data cards so that successful operation of the program will be achieved. 


DATA SET 

The data set consist of 11 input parameters that must be properly presented on the 
data cards in order to make proper use of the program. The parameters are defined 
as follows: 

1. R = radius of curvature in mm. 

2. K = order of diffraction, where K = ±1,±2, etc. The choice of pegative or posi- 
tive orders is at the discretion of the program user. 

3. N = groove frequency in grooves/mm. 

4. X Q = laser recording wavelength in angstroms. 

5. X* = correcting wavelength in angstroms. 

6. Xl = lower wavelength limit in angstroms for displaying the performance of the 
grating. 

7. AX = wavelength interval in angstroms for displaying the performance of the 
grating. 

8. Xu = upper wavelength limit in angstroms for displaying the performance of the 
grating. 

9. a L = lower angle of incidence in degrees decimal form. 

10. «u = upper angle of incidence in degrees decimal form. 

11. A a = incremented angle of incidence in degrees decimal form. 

a L> a u> an( i Aa require some explanation. The incremental angle of incidence Aa 
provides the program user the opportunity to obtain optimum recording parameters for 
successive angles of incidence at specified increments without punching data card sets 
for each desired angle of incidence. For example; if a L a °d “U are g» ven respective 
input values 0.0° and 10.0° and Aa = 5.0°, then optimum recording parameters will be 
obtained for a = 0.0°, 5.0°, and 10.0°. Complete data will be displayed for each of 
the three angles of incidence. If the program user selects one input angle of incidence, 
then ay and ay l' ave *h e 53016 input va l ue a °d Aa = 0.0°. 
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DATA CARDS 


The order of the first eight input parameters must be punched in the order shown in the 
DATA SET section starting in column 1 in format F7.1 with seven columns reserved for 
each of the eleven parameters on the first data card and each parameter must include a 
decimal. The parameters a U> an ^ are punched on the second data card in that 
order with seven columns reserved for each of these three parameters in F7.1 format and 
each parameter must include a decimal. Figure 3 shows a typical two card data set to 
illustrate the description stated above. 


/ 1000.0 TO 1200 0 4070 3 2000.0 1000 0 000 0 4000.0 

R K N ,\ 0 X* X L \\ \ 0 

to o g g g g g g r iu g o o o g g o j a o i g g o g g g g i o u g g o c a g g g o o g # c g g g o o g i 
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Figure 3. 


NOTE: 


1. The program user may include as many data sets as required in a single run 
of the program. 

2. The values of all 1 1 input parameters must be punched on all data card sets. 
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PROGRAM EXAMPLE AND DATA DISPLAY 


This section presents a typical example of the use of the program and also describes the out- 
put data in detail. 

Example: Determine the optimum recording parameters 7 and S such that astigmatism or 
coma is made zero at 2500A for angles of incidence a = 0.0° to 25.0° in the incremental 
angles of 2.0°. The grating has 1200 L/mm with a radius of curvature of 1000 mm and will 
operate in the 1 st order. The wavelength performance range is from 1500A to 4000A in 
500A intervals. 

Solution: A data card set is punched using the data supplied in the example and punched 
in the order listed as follows: 

First Card 


1. R= 1000.0 mm 

2. K = +1 .0 

3. N = 1200 grooves/mm 

4. X 0 = 4579.3A 

5. X* = 2 5 00.0 A 

6 . X L = 15 00.0 A 

7. AX = 500.0A 

8 . Xu = 4000.0A 
Second Card 

9. a L = 0.0° 

10. <*u = 25.0° 

11. Aa = 2.0° 

The output data for the example cited is shown in Figure 4. Output data for only one 
of the angles of incidence will be presented for the purpose of illustration. 

A detailed description of the output data for correcting astigmatism or coma at 2500A 
is as follows: 
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Figure 4. (Continued) 






Line 1 is self explanatory and refers to the next two lines. 

Line 2 shows the equations for calculating L^, Lg, Lq and Lq. 

Line 3 shows the radius of curvature R, recording wavelength X 0 = RLAM, grating order 
K, and the groove frequency N. 

Line 4 shows that astigmatism is being corrected and is indicated by option = 1. Option 
= 1 is an internal part of the program and is not an input parameter. 

Line 5 shows the number of solutions available at the correcting wavelength X*. 

Line 6 indicates that the following data refers to the first solution. 

Line 7 shows the values of the correcting wavelength X* = LAMDA, angle of incidence 
a = ALFA, diffracted image angle 0 = BETA, recording angle y = GAMA and re- 
cording angle 6 = DELTA. 

Line 8 shows the values of the aberrations corrected at X*. Note the value for astigma- 
tism A2 is essentially zero. 

Line 9 shows the values of the aberrations for the conventional equivalent grating at X*. 
Notice how much greater the value of astigmatism is at X* compared to the value 
for the holographic grating in Line 8. 

Line 10 indicates that the data in Lines 11 and 12 refer to the holographic grating. 

Line 1 1 shows the titles for the performance of the holographic grating. LAMDA is the 
wavelength at which the grating is evaluated relative to the aberrations. BETA is 
the diffraction angle. The aberration titles Aj , A2, Cj , C2, Sj , S2 and S3 have 
been defined in the text. 

Line 1 2 shows the values of the wavelength, diffraction angle, and aberrations for the 

performance of the grating. Note the values of the aberrations at 2500A and com- 
pare with Line 8. 

Line 1 3 indicates that the table that follows refers to the performance of an equivalent 
conventional grating. 

Lines 14 and 15 refer to the conventional grating. Compare the values of this table with 
the values in the holographic grating performance table. Note the reduction in 
the aberration values for the holographic grating. 

Lines 1 6 indicates that the following data is for correcting coma and is indicated by op- 
tion = 2. Option = 2 is inherent to the program and is not an input parameter. 
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Line 17 indicates that for the angle of incidence shown in Line 19, there are two solu- 
tions for correcting coma. 

Line 18 shows that the first solution will be evaluated for coma. The output data for 
the second solution is not shown. 

Lines 19 through 27 are defined in the same way as Lines 6 through 15 except that they 
refer to coma. 
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SIMULTANEOUS CORRECTION OF ASTIGMATISM, COMA, AND 
FIRST SPHERICAL ABERRATION TERM OF HOLOGRAPHIC 
CONCAVE GRATINGS FOR ROWLAND CIRCLE 
SPECTROGRAPHS 


SIMULTANEOUS CORRECTION OF ASTIGMATISM, COMA, AND 
FIRST SPHERICAL ABERRATION TERM OF HOLOGRAPHIC CONCAVE 
GRATINGS FOR A ROWLAND CIRCLE SPECTROGRAPH 


INTRODUCTION 


When designing spectroscopic instruments such as spectrographs, monochromators, spec- 
trophotometers and other instruments employing concave optical reflecting diffraction 
gratings, one of the most important optical design criteria is to design the optical system 
such that optical aberrations are reduced or eliminated over the desired wavelength range 
of operation. The most notorious aberration is normally astigmatism because the dif- 
fracted images are magnified in the vertical direction. As a result, large area detectors are 
required to sense all the energy in the diffracted image or, as is the usual case, available 
small area detectors are used at the expense of losing a great deal of energy in the dif- 
fracted image. 

However, with the advent of aberration corrected concave holographic optical reflecting 
gratings, astigmatism and other aberrations have been reduced or eliminated at specified 
wavelengths. Recent publications have shown that for certain spectroscopic instrument 
designs incorporating corrective type holographic gratings either astigmatism or coma may 
be minimized over a desired wavelength range or astigmatism or coma may be eliminated 
at specified wavelengths/ 1 . 2 . 3 , 4 ) However, the aberrations can not be corrected simulta- 
neously at specified wavelengths. The purpose of this paper is to present the theory and 
basic equations that lead to the design of concave holographic diffraction gratings for 
Rowland circle spectrographs that do eliminate astigmatism, coma, and the first spherical 
aberration term simultaneously at specified wavelengths. However, the simultaneous eli- 
mination or correction to zero of these aberrations occurs at wavelengths that are greater 
than or equal to the laser recording wavelength. 

A corpputer program was also developed to numerically analyze the theory, determine 
the optimum instrument performance and recording parameters necessary for grating 
fabrication. The program also displays the performance, for comparison, of both the 
holographic grating and equivalent conventional grating as a function of wavelength. 


THEORY AND BASIC EQUATIONS FOR DESIGNING HOLOGRAPHIC SPHERICAL 
CONCAVE GRATINGS 


Figure 1 shows the optical system for using and recording a concave spherical holographic 
grating. The rectangular coordinate system is defined as follows: I is the origin at the 
vertex of the grating blank and the x-axis is normal to I. The x-y plane is defined l?y I 
and the two recording point sources C and D. Points A, M, and B are, respectively, a 
self luminous point on the entrance slit, a point on the grating surface, and a point on 
the diffracted image from M at a wavelength X in the n th order. A knowledge of the 
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Figure 1 . Schematic diagram of the optical system. A is a self luminous point on the 
entrance slit, M is a point on the n th groove and B is the spectral image point. C and 
D are the recording point light source positions. The performance and recording param- 
eters are defined as follows: 

Performance Parameters : 

L a is the distance from A (entrance slit) to the center of I of the grating. 

Lg is the distance from the diffracted image B to the center I of the grating. 
a is the angle of incidence. 
p is the diffraction angle, 
z is the object height. 
z is the image height . 

Recording Parameters: 

Lc is the distance from the point source C to the center I of the grating. 

Lq is the distance from the point source D to the center I of the grating. 

7 is the position angle for point source C relative to the normal to the grating. 

6 is the position angle for point source D relative to the normal to the grating. 
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geometry of the optical system, Figure 1, leads to writing the aberrant optical path func- 
tion AM for the conventional grating as 

AM = AM + MB - (IA + IB) - KNX (1) 

where K is the order of diffraction, X is the diffracted wavelength and N is the grating 
groove frequency. 

When a grating blank coated with a photoresist is illuminated by coherent spherical wave 
fronts originating from points C and D, grooves are recorded, after development, in the 
photoresist in response to the interference fringes produced at the site of the grating 
blank. The interference fringes are formed according to the equation 

N -! ' (2) 


where X 0 is the recording laser wavelength. 

Substituting eq. (2) into eq. (1 ) gives the aberrant optical path function 

AM = MA + MB-(IA + 1B)-^ j(MC - MD) - (IC - ID) J (3) 

lor the holographic grating. 

The quantities on the right side of AM in equation (3) are derived based on the coordinate 
system described in Figure 1 and the equation of a sphere. Since these derivations are 
lengthy, we suggest the reader see reference 4 for a more detailed analysis. Therefore it 
suffices to say that after a series expansion is performed, the following result for AM is 
obtained. 


AM = -Y[sina + sin/3-^ | sin y -sin 5 } ]+~[Aj] (4) 

+X l A 2l + T ICll + -T^IC 2 ) +Y 4 [Sj] 

+ Y 2 Z 2 [S 2 ] + Z 4 [S 3 ] + ,.. 

where Y is the width and Z is the height of the grating and 

a - c° s2 <* j. cos 2 j3 cos /3 KX [ cos 2 7 cos-y / cos 2 5 cos 6 \ | 

'' L A --r + l„ - r (tj — ir*l-T5 — r-/) 

•a _ I cosfl , 1 cos/3 KX f 1 cos-v / 1 cos6\l 

A2 ~ it ' tt + 1; - R -T£{q:--r*(r5' “)) (5) 
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Note, that in the expressions for equation ( 5 ), L A , Lq, a, and j 3 are the terms for the 
conventionally ruled grating and Lc, Lq, 7, and 6 are the terms related to the hologra- 
phic contribution. R is the radius of curvature of the grating and K is the order. A 0 js 
the recording laser wavelength. The expressions in equation ( 5 ) are defined as follows: 
Aj is the horizontal focusing term, A2 is the astigmatism or vertical focusing term, Cj 
and C 2 are coma terms and Sj , S2, S3 are spherical aberration terms. 

The grating equation is given by 


a (sin a + sin 0) = KX 

and for the holographic grating, 0 is the effective grating constant defined by 

a ~ 01 - c where 7 > 6 
sin 7 - sin 0 


( 6 ) 


( 7 ) 


The above description presents the general theory and basic equations for designing aber- 
ration corrected holographic gratings. We will now apply the theory to design a holo- 
graphic grating that simultaneously corrects astigmatism, coma, and the first spherical 
aberration term at a specified wavelength. 


HOLOGRAPHIC GRATING DESIGN FOR ROWLAND CIRCLE SPECTROGRAPH 


Aj and Cj in equation ( 5 ) must be made zero in order to design a holographic grating 
that has the Rowland circle as its focal curve. The solution for this condition is: 

L A = R cos a, Lg = R cos ft L(; = R cos 7 and Lq = R cos 5 (8) 

Substitution of equation (8) into the aberrant optical path function reduces equation ( 4 ) 
to 

AM = -Y[sin a + sin fi - ^ | sin 

+ 2rT IC 2 ) 1 + W S 2 1 + 8R? S 3 1 


7 -sin 6 ] +^(A 2 '1 


where the expressions in equation ( 5 ) take the form 



sin2 q* 
cos a* 


+ 


sj n 2 ft* KX* 
cos/J* ~ Xq 


sin 2 7 _ sin 2 8 
cos 7 ~ cos 5 


( 10 ) 
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c 2 ' = 

sin 3 a* 

sin 3 0 * 

KX* | 

| sin 3 7 sin 3 5 

cos 2 a* 

cos 2 0 * 

Xo i 

| cos 2 7 cos 2 5 

V - 

sin 2 a* 
cos a* 

sin 2 0 * 

+ - 

COS0* 

KX* j 
Xo 1 

1 sin 2 7 sin 2 6 | 

| cos 7 cos 8 ) 

sv = 

sin 2 a* 

(1 + sin 2 a*) 

sin 2 

4. -- 

0 * (1 + sin 2 0 *) 

cos 3 a* 


cos 3 0 * 


KX* j 

sin 2 7(1+ sin 2 7) 

sin 2 6 ( 1 + sin 2 6 ) 


" x G I 

cos 3 7 


cos 3 5 
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sin 2 a 

-s „ . . sin 2 0 * 
(2 cos 2 a* - 1 ) + ~ 

cos J a 

cos 3 0 * 

KX* ( 

sin 2 7 - 

- — ^ 

— (2 cos 2 7- 1) - 

I 

cos 3 7 


sin 


cos- 


i£ 
3 6 


(2 cos 2 5 - 1 ) | 


X* is the wavelength at which the aberrations are corrected to zero and 0* is calculated from 


a (sin a* + sin 0 *) = KX* 

Note that the groove frequency of the cratine is given by N = — . 

a 


( 11 ) 


SIMULTANEOUS CORRECTION OF ASTIGMATISM, COMA, AND FIRST SPHERICAL 
ABERRATION TERM AT X* 


The simultaneous elimination of astigmatism, coma, and first spherical aberration term at 
X* depends on the instrument and recording parameters. For given instrumental param- 
eters o and X* there are specific angles of incidence a* for which the recording angles y 
and 6 are chosen such that astigmatism, coma, and first spherical aberration term are 
simultaneously made zero. These aberration corrections are accomplished by simultane- 
ously solving the following set of equations: 

a (sin a* + sin 0*) = KX*, ( 1 1 ) 

- X 0 

0 ~ sin 7 - sin 5 (7) 
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and A 2 ' = C 2 ' = 0 in equation (10). Since the expression fqr Sj ' is identical to A 2 *, the first 
spherical aberration term Sj is automatically zero. A numerical analysis of the above theory 
shows that solutions exist only for X*> X 0 . The numerical analysis also shows that there 
are one or more angles of incidence a* for each correcting wavelength X*, depending upon 
groove frequency for which A 2 ' = C 2 ~ V = 0. Figure 2 shows some of a family of 
curves which for a given groove frequency and correcting X*, the angles of incidence a* 
and the recording angle 6 may be obtained such that A .2 ~ = Sj' = 0. The groove 

frequency is shown for each curve and the correcting wavelengths are indicated at differ- 
ent points along that curve. For example, at a. groove frequency of 1200 L/mm and a 
correcting wavelength at 5000 A, there are three angles of incidence (5.07°, 30.77° and 
43.16°) and two corresponding recording angles 6(-2.83° and 7.01°) where the aberrations 
may be simultaneously corrected. It is interesting to note that as the grating groove fre- 
quency increases, the number of incident angles a* and range of correcting wavelength 
decreases. 

The diffracted angle 0* and other recording angle y are calculated by equations (1 1) and 
(7) respectively. After the instrument and recording angles have been determined, the 
remaining instrument parameter values La, Lg, and recording parameter values Lc, Lp 
are calculated using expressions in equation (8) along with the value of the radius curva- 
ture R. 


COMPUTER PROGRAM 


In the first part of this paper, we have presented the theory and basic equations necessary 
to design holographic gratings that are aberration corrected at specified wavelengths at or 
above the laser recording wavelength for Rowland circle spectrographip mounts. We have 
also presented data (Figure 2) that shows the different angles of incidence that could be 
used to design an optical system based upon a specified groove frequency and correcting 
wavelength that eliminates astigmatism, coma, and first spherical aberration term simul- 
taneously. Now we present the computer program that made it possible to numerically 
analyze the theory and calculate the optimum performance and grating fabrication param- 
eters. The program also displays the performance of the holographic grating and equiva- 
lent conventional grating relative to all the aberrations as a function of wavelength. 


DESCRIPTION OF COMPUTER PROGRAM 


The computer program consists of a main program and four subroutines. The main pro- 
gram directs the flow of calculations and calls upon the subroutines to perform their 
various functions. The flow diagram (Figure 3) shows the events that will be performed. 

The algorithm may be summarized as follows: Assume that the radius of curvature, groove 
frequency, and correcting wavelength has been selected, then an angle of incidence oq * 
0.0° is chosen and the corresponding diffracted angle 0 is calculated. Based on the above 


7 


n oo? i 



Figure 2. Family of curves showing the wavelengths in nanometers at which A 2 = C 2 = Sj = 0. The correcting wavelength in nm 
is along each curve and is a function of the recording angle 6 and the incident angle a. The groove frequency is shown at the end 
of each curve. 





Figure 3. Program Flow Diagram. 
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parameters, the recording parameters and 8j are determined such that A2 = 0 . The 
above procedure is repeated at another angle of incidence 0:3, which is determined by an 
equation of the form 0:3 = oq + Aa; where aj is the previous angle of incidence and Aa 
is a very small incremental angle. At this point, there are two sets of parameters, oq , 0 \, 
7 j, 6j , and aj, 02 , To, 63 f° r which A2 = 0 . Note that A2 can be made zero at any 
angle of incidence. For these two sets of parameters, two coma values for C2 correspond- 
ing to Coloq ) and €2(0:2) are calculated. If the values of the C2 terms have opposite signs, 
then there exists a value a between Ofj and 03 for which A2 = C2 = 0 for the same values 
of a, 0 , 7, and 6. The values of a and associated parameters are determined by a modifi- 
cation of Newton’s method. 

However, if C2(<*i) and 02(0:2) have the same sign then there isn’t any value of a between 
oq and a.2 for which A2 = C2 = 0 . When this happens, we set 0:3 = 0:2 + Aa and calculate 
the corresponding angles /J3, 73, and 83 such that A2(at3) = 0 . Then 03(0:3) is calculated 
and the sign is compared with 03(0:2). Again, if the signs are opposite there is another 
value of a between <*3 an< * a 3 such that A2 = C2 = 0 . This process is repeated until all 
angles of incidence between 0 ° and 90 ° have been evaluated and only those angles of in- 
cidence are selected which satisfy the condition that for a certain 0; there is a 0 , 7, and & 
that produces A2 = C2 = 0 for a specified groove frequency, radius of curvature and cor- 
recting wavelength. Since the equations for A2 and Sj are identical, then Sj is also zero. 

The algorithm is supported by four subroutines. Subroutine ASTGSM calculates the 
parameters required to make astigmatism zero at angles of incidence a. 

Subroutine FUNCAC calculates the values of the coma terms based on the parameters 
used to make A2 = 0 . Then it compares the signs of the coma terms; if opposite this 
subroutine will determine the angle of incidence and related parameters that make coma 
zero. If the signs are the same, it returns the problem to subroutine ASTGSM and the 
process is repeated. 

Subroutine Coma calculates the values of a, 0 , 7 and 8 such that C2 = 0 , when A2 - 0 
for the same set of parameters. 

Subroutine PRFORM calculates and displays in the form of a table the numerical values 
of each aberration for the holographic grating and equivalent conventional grating as a 
function of wavelength. 


INPUT DATA, DEFINITION OF TERMS AND DATA CARDS 


This section describes the input data and defines the parameters that will be punched on 
the data cards so that successful operation of the program will be achieved. 


(A) DATA SET 

The data set consists of 1 2 parameters that must be properly presented on the data cards in 
order to successfully use the computer program. These parameters are defined as follows: 
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1. R = radius of curvature in mm. 

2. N = groove frequency in grooves per mm. 

3. X 0 = laser recording wavelength in angstroms. 

4. X* = correcting wavelength in angstroms. 

5. Xl = lower wavelength limit in angstroms for displaying the performance of the 
grating. 

6. Xjj = upper wavelength limit in angstroms for displaying the performance of the 
grating. 

7. AX = wavelength interval in angstroms for displaying the performance of the grating. 

8. K = order of diffraction, where K = ±1, ±2, etc. The choice of negative or posi- 
tive orders is at the discretion of the program user. 

9. Y = grating width in mm. 

10. Z =, grating height in mm. 

11. Aa = incremental angle in degrees decimal form. The smaller the incremental angle, 
the more accurate the successive values of a. 

12. lOPTION = Selects whether or not the coefficients Y and Z should be included in 
the aberration values. If lOPTION = 1 is selected, then the aberrations are calcu- 
lated without the coefficients. If lOPTION = 2 is selected, then the aberrations 
calculated include the coefficients. 


(B) DESCRIPTION OF DATA CARDS AND SAMPLE DATA CARDS 

The order of the first eight input parameters must be punched in the order shown in 
(A) starting in column 1 in FI 0.1 format with 10 columns reserved for each of the 
seven parameters and each parameter must include a decimal. The eighth parameter K 
is in 12 format and punched in columns 71, 72. The next four parameters are punched 
on the second card where parameters Y, Z and Aa are format FI 0.1 and lOPTION is 
in 12 format. Figure 4 shows a typical two card data set to illustrate the description 
stated above. 
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• Figure 4. 

NOTE: 1. The program user may include as many data sets as required in a single 

run of the program. 


2. The values of all 12 input parameters must be punched on all data card 
sets. 
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EXAMPLE AND DESCRIPTION OF DATA 


A typical example of the program and .detailed description of the output data is as follows: 

Example: Design a holographic diffraction grating that simultaneously corrects for astig- 
matism, coma, and the first spherical aberration term at 5000 A for a Rowland circle 
spectrograph. The grating has 1200 L/mm and a 1.0 meter radius of curvature. Deter- 
mine the performance of the holographic grating and equivalent conventional grating 
throughout the wavelength range 4500 A to 9000 A at 500 A intervals. 

Solution: Data cards are punched using the data supplied in the example and punched 
in the order listed below: 

First Card 


1. 

R = 1 000 mm 

2. 

N = 1200 L/mm 

3. 

X 0 = 4579.3 A 

4. 

A c = 5000 A 

5. 

X y = 9000 A 

6. 

X L = 4500 A 

7. 

AX = 500 A 

8. 

K = +1.0 

Second Card 

9. 

Y = 1.0 

10. 

N 

It 

o 

11. 

o 

to 

o 

II 

< 

12. 

IOPTION = 1 


The output data for the example cited is shown in Figure 5. Output data for only one 
of the angles of incidence will be presented for the purpose of illustration, 

A detailed description of the output data for simultaneously correcting astigmatism, coma, 
and the first spherical aberration term at 5000 A is as follows: 
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Line 1 : Shows the correcting wavelength LAMDA = X c , angle of incidence ALFA = a, 

diffracted wavelength BETA = j3, recording angle GAMA = 7 and recording 
angle DELTA = 5. 

Line 2: Self-explanatory. 

Line 3: Shows the symbols for the wavelength and each aberration. The symbols for 

the aberrations have been defined in the text. 

Line 4: Shows the correcting wavelength and values for each aberration at the correcting 

, wavelength. 

Line 5: Self-explanatory. 

Line 6: Indicates that the holographic grating has been evaluated. 

Line 7: Shows the titles of wavelength = LAMDA and the aberrations Aj, A2, Cj, C2, 

Sj, S-> and S3. The table shows the values of each aberration as a function of 
wavelength. 

Line 8: Indicates that Line 9 and the corresponding table are for the evaluation of the 

equivalent conventional grating. The equivalent conventional grating is defined 
as having the same groove frequency and radius of curvature as the holographic 
grating. However, the conventional grating cannot correct aberrations. Note 
the reduction in aberration values as a function of wavelength as compared to 
the conventional grating. For example, at 5000 A the value of astigmatism A2 
for the holographic grating is 0.316973 x 10"4 and that for the conventional 
grating is 0.31292. This shows that astigmatism is significantly reduced by a 
factor of ! Q 4 for the holographic grating at 5000 A. Note that C2 and Sj are 
also reduced by the same order of magnitude. The remaining spherical aberra- 
tion terms are also reduced by a factor of 10 2 for the same wavelength. Also 
note that the aberrations increase for decreasing and increasing wavelengths 
about the correcting wavelength at 5000 A. However, the increase in aberra- 
tions remains less than that of the equivalent conventional grating. 

In summary, the theory and basic equations have been presented to design concave spher- 
ical holographic reflecting diffraction gratings that are corrected to zero at specified wave- 
lengths at or above the laser recording wavelength. A numerical analysis shows that there 
are specific angles of incidence based on the selected groove frequency and correcting 
wavelength that simultaneously corrects the aberrations A2, C2 and Sj. A computer pro- 
gram lias been developed so that application of the theory may be applied to practical 
spectroscopic instrument design requirements for space and ground based laboratory 
applications. 
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DESIGN OF ABERRATION CORRECTED HOLOGRAPHIC TOROID 
GRATINGS FOR SEYA-NAMIOKA TYPE MONOCHROMATORS 


Part I: Minimization of Aberrations For 
Desired Wavelength Range 

A. J. Caruso 

M.S. Bhatia 

Seya-Namioka type monochromators are designed such that wavelength selection is obtained by 
rotating a concave reflecting grating about its vertical axis. Diffracted image focussing is achieved 
by one optical element; normally a concave spherical reflecting diffraction grating. The angle of 
deviation, defined as the angle between the entrance and exit slits relative to the vertex of the 
grating, is a fixed value, normally 70°. At this fixed value most of the optical aberrations are 
reasonably acceptable except for severe astigmatism which is caused by the mechanically ruled 
concave spherical grating. With the advent of ruling holographic aberration corrected concave 
spherical reflecting diffraction gratings, it has become possible to design Seya-Namioka type mono- 
chromators such that astigmatism or coma may be corrected to zero at specified wavelengths or 
minirnized over a desired wavelength range (ref, 1 , 2 , 3 , 6 ). These holographic grating designs 
have resulted in a significant improvement for the output efficiency of Seya-Namioka monochro- 
mators for wavelengths longer than 100 nm. However, for spherical concave gratings with angles 
of deviation greater than 70° holographic grating theory shows that aberration correction cannot 
be achieved for wavelengths shorter than 80 nm. Further, at 70° angle of deviation the output 
efficiency is also decreased because grating coating materials reflectivity is significantly reduced 
at wavelengths below 10.5 nm. 

In order to over come this problem, we have substituted the toroid geometry for the grating 
blank in place of the concave spherical geometry, (reference 5) The major advantage of the toroid 
geometry is that astigmatism can be eliminated at specified wavelengths and also dramatically mini- 
mized over a wide wavelength range by properly determining the correct toroid vertical radius of 
curvature and horizontal radius of curvature as a function of the angle of incidence and the required 
holographic ruling geometry. Another advantage is that at larger incident angles the reflectivity 
of the same grating coating material increases at wavelengths less than 10.5 nm. 

The purpose of this report is to present the theory and design equations for holographically ruled 
toroidal diffraction gratings that minimize aberrations for Seya-Namioka type monochromators 
with angles of deviation larger than or equal to the normal 70° for wavelengths shorter than 100 
nm. 


This report is divided into four sections. Section 1 presents the theory and design equations for 
designing aberration corrected holographic toroidal gratings. Two methods have been developed: 
the general method requires that the instrumental parameters be related only to the holographic 
toroid grating design. The second method is called the modified method. This method requires . 
that the instrumental parameters be the same for both the conventional and holographic toroidal 
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gratings. The modified method therefore designs such that the holographically ruled toroidal 
grating is interchangeable with the conventionally ruled grating providing they have the same 
groove frequency, and radius of curvatures R and p t . 

Section 2 presents a description of the computer program. The program has been developed such 
that for a given groove frequency and wavelength range, the optimum instrumental and holographic 
ruling parameters are determined for the toroidal concave grating. The program also displays the 
performance of the holographic toroidal grating and equivalent conventional toroidal grating, as 
a function of wavelength relative to optical aberrations such as astigmatism, coma, and spherical 
aberrations. Section 3 presents the input data, definition of terms and data card description. 
Section 4 presents an example for using the program and describes in detail the output data. 

Section 1 

Theory and Basic Equations for Designing Holographic Toroidal Concave Gratings 


A toroid is produced when a circle is rotated about an axis which lies in the plane of the circle but 
does not pass through the midpoint of the circle (Fig. 1). Figure (1) also shows a section (toroidal 
grating blank) cut out of the toroid for which there are two radii of curvatures: a long horizontal 



Figure 1. Schematic diagram of a toroid. R and p t are the radii of the toroidal surface. Cut 
out section is typical of a toroid grating blank 
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radius R and a short vertical radius p t . If p t = R, then a spherical concave surface will be gener- 
ated; ie, when the axis of rotation passes through the midpoint of the circle. As a result, all infer- 
ences which are valid for the toroid are also valid for the sphere when p t = R. 

We now proceed with the theory and basic equations for designing holographic toroidal gratings 
with the above concept in mind. 

Figure (2) shows a toroidal grating blank in an optical schematic diagram which defines a rectan- 
gular coordinate system. Let the origin be at vertex 0 of the concave toroid grating blank which 
has a horizontal radius of curvature R and a vertical radius of curvature p,. Let the x-axis be de- 
fined by O and two coherent point light sources C and D which are used to produce the inter- 
ference fringes on the concave toroidal blank surface. Points A, P, and B are self luminus points 
on the entrance slit, a point on the grating and a point at the focused diffracted image from P at 
wavelength X in the order respectively. A knowledge of the geometrical optical system of 
Figure 2 leads to writing the aberrant optical light path function F for the conventional grating 
as 


F = AP + PB + nin X ( 1 ) 

where n is the groove frequency, m is the diffracted order and X is the diffracted wavelength. 

When a grating substrate coated with a photoresist layer is illuminated by coherent spherical wave- 
fronts originating from points C and D, grooves are recorded and formed after development in 
the photoresist in response to the interference fringes produced at the site of the grating surface. 
The interference fringes are formed according to equation 


nX 0 = [(CP) -(DP)] - [(CO) -(DO)] 
where n is the groove frequency and X 0 is the laser wavelength radiation. 


( 2 ) 


Substituting equation (2) into equation (1) gives the aberrant optical light path function 


rroA 

F = AP + PB +—- [(CP - DP) - (CO - DO)] 

Xn 


(3) 


for the holographic grating. The terms on the right side of F in equation (3) are derived based on 
the coordinate system described in Figure (2). By substituting the trignometric terms for the 
toroid and other required terms, in equation (3), the following aberrant optical light path function 
F is derived after a power series expansion is performed. 


F = F 000 +WF 100 + LF 011 +‘/ 2 W 2 F JOO +V 2 L 2 F 020 
+ , /iW 3 F 300 +'/:WL j F I20 +WLF,,, 

+ 1/8 W 4 F 400 + !4 W 2 L 2 F 220 + 1/8 L 4 F 040 '/4WL 2 F 202 
+ */4L 2 F 022 +'/2L 3 F 0 3, +V4W 2 LF ail +.... 


(4) 
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Figure 2. Schematic diagram showing the geometry for recording the holographic grating. The 
instrumental and recording parameters are defined as follows: 


Instrument Parameters 

r 0 is the object distance from the entrance 
slit to the center of the grating. 

rjj is the image distance from the exit slit to 
the center of the grating. 

a is the angle of incident radiation relative 
to the normal of the grating. 

0 is the diffracted image angle relative to 
the normal of the grating. 

A is the entrance slit position and B is the 
exit slit position. 


Recording Parameters 

is the distance of recording point source C to 
the center of the grating. 

tp is the distance of the recording point source D 
to the center of the grating. 

6 is the recording angle for positioning point 
source D relative to die normal of the grating. 

7 is the recording angle for positioning point 
source C relative to the normal of the grating. 
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(5) 


♦ and 

1 ^ijk “ Mijk +^~ Hjj k 

where the subscripts ijk of Fjj k are exponents ofW'UZk except for F 0 co- 
The terms M ijk and Hjj k are defined in terms of the function fjj k as follows 

M ijk = f.jk (P. °0 + t'ijk (p'. 0) 

Hjjk “ Ojk (Pc, 7) — fjjk (Pd » fi) 

for ijk = (200) and (300) 
and 

Mijk = fjjk (P. Q. <*) + fijk (p'. q, 0) 

Hijk = ^ijk (Pc » q» 7) Ojk (Pd > Q> 

for ijk = (020) and (120) 

M ijk = M 0 J0 -(-!)** [fijk (P. <*) + fijk (p',0)l 
Hijk =H 020 -C-l) 444 If Uk (Pc, 7) "fijk (PD,*)] 

for ijk = (400) 
and 

Mjj k = M 020 -(-l)* 4 [f ijk (p, q,«) + f ijk (p', q,0)l 
HUk = M 02 o-(-l) i4i [fijk (Pc, q, 7) "fijk (p D ,q,fi)] 
for ijk = (220) and (040) 

where the expressions for the conventional toroid grating are 

fj 0 0 (P, °0 + f(P*, 0) = cosa (pcosa-1) + cos/3 (p'cos0-l) 

fo 2 0 (P, q, a) + f (P'> q, 0) = P “ qcosa + p' - q cos/3 

fj 0 0 (P» a) + f (p', 0) = psinacosa (pcosa - 1) + p'sin/3cos/3 (p'cos0~l) 

f, jo (p, q, a) + f (p' q, 0) = psina(p-qcosa) + p' sin/3 (p'-qcos0) 

f 40 o (p, cr) + f 40 o (p',0) = p - cosa-pcosa (pcosa-1) [cosa(pcosa-l) - 4 psin 2 a] 


(6a) 

(7) 


(7a) 


(8a) 

(8b) 

(8c) 

(8d) 

(8e) 
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+ p' — cos /3 - p' cos /3 (p' cos/ 3 - 1 ) [cos /3 (p' cos/ 3 - 1 ) — 4 p'sin 2 / 3 ] 
f 22 o (P. q. «) + ^220 (P> q, / 3 ) = q (psin 2 o!- cosa) + p 2 cosa 

+ ( 3 p 2 sin 2 a-p 2 )(p-q cosa) (80 

+ q (p'sin 2 /3 — cos/ 3 ) + p' 2 cos /3 + ( 3 p' 2 sin 2 / 3 -p' 2 ) (p'-q cos/ 3 ) 


l*o 4 o (P. q.«) + f*o4o (P'. q. 0) = q 2 (p-qcosa) - p(p-qcosa) 2 

(8g) 

+ q 2 (p' - q cos/ 3 ) -p' (p'-q cos/ 3 ) 2 
and the expressions for the holographic contribution are: 

t'200 (Pc> 7 ) - 1*200 (Pd . 5 ) = cost(p c cos t“ 1) - cos 5 (p D cos 5 -l) (8I1) 

1*020 (Pc. q, 7 ) “ t'020 (Pd. q. 5 ) = P C “ qcos7~p D -q cos 5 (8i) 

Ooo (Pc. 7 ) - 1*300 (Pd > 7 ) = Pc sin 7 cos 7 (p c cos7-l)-p D sin 5 cos 5 (p c cos 5 - 1 ) (8j) 

1*120 (Pc, q, 7 ) ~ O20 (Pc, q. 8) = P c sin 7 (p c -q COS7) - [p D sin 5 (p D -q cos 6)] (8k) 


1*4 0 0 (Pc, 7) ~ 1*400 (Pd. 5) = Pc cos 7 — p c cos7 (p c COS7- 1) 

(cos 7 (p c cos 7 - 1) - 4 p c sin 2 7) 

(81) 

- [Pd “ cos 5 -p D cos 5 (p c cos 7 ~ 1) 

(cos 6 (p c cos 5 - 1) - 4 p c sin 2 5) 

1*2 20 (Pc.q.7)-f2 20 (Pn.q. 8) * q (Pc sin 2 7 — cos 7) + p 2 cos 7 

+ (3p 2 sin 2 7 ~ Pc ) (Pc - q cos 7) (8m) 


“ [q(PD sin 2 5 - cos 5) + p£ cos 6 +(3p£, sin 2 5-p 2 )(p D -q cos 6)] 
1*040 (Pc, q. 7 )~ f*040 (Pd, q. 8) = q 2 (p c ~ q cos 7) - p c (p c -q COS7) 2 


~ [q 2 (Pd ~ q cos 8) -p D (p D -qcos 5 ) 2 I 


(8n) 


where 

R , R R R 

P — > P — 1 1 Pc ~ » Pd — — 

r r r c ro 


( 9 ) 


6 


and q= — where R is the horizontal radius and p. is the vertical radius of the toroid. 

Pt 

The term F, 00 in equation 4 is the well known grating equation and is given as 

mX = a (sin a + sin 0) ( 1 0 ) 

for the principal ray, the ray originating from the center of the entrance slit and diffracted from 
O at an angle 0 . In order to holographically rule the grating with the required groove frequency, 
the effective grating constant is defined as 


sin 5 - sin 7 

where 5 > 7 . 7 and 5 are the recording angles. 

Note that the number of grooves per unit length is n =—. 

o 

The individual expressions for equations ( 8 a) thru ( 8 n) have physical significance as to the forma- 
tion and imperfections of the diffracted image. The expressions control the various conditions 
for the foci and image deficiencies in the following order; 8 a and 8 h govern the primary focusing 
condition on the Rowland Circle, 8 b and 8 i determine the amount of astigmatism relative to the 
diffracted image, 8 c, 8 j, 8 d, and 8 k Coma, 8 e and 81 higher order focusing condition on the Row- 
land Circle and 8 f, 8 g, 8 m and 8 n the tone aberrations. It is interesting to note that when one 
compares these expressions of the toroid with those for the equation of a sphere Ref (1), there 
are some expressions that are identical for both, (toroid & sphere). For example, all the expressions 
which do not include q are identical with those expressions for the sphere. It is only those ex- 
pressions that involve q that account for the difference between the two gratings. Note that the 
difference between the two gratings vanishes when p t = R. 

Minimization of Aberrations 


The theory described up to this point has been presented for the general case. We now apply the 
general theory to a specific application: The Seya-Namioka monochromator. This type of mono- 
chromator has a single reflecting and wavelength dispersing element; the grating. Wavelength selec- 
tion is achieved by rotating the grating about the vertical axis at the front surface of the center of 
the grating. The basic instrumental parameters for this type monochromator are defined as follows; 
see fig 3. 


r 0 , fo ~ constant, 2 K -a 0 -p o - constant 
a 0 = K + 0, 0 O = 0 - K 


( 12 ) 


where r 0 is the distance from the center of the extrance slit to O and rj, is the distance from O to 
the center of the exit slit. 2K is the angle of deviation AOB and 0 is the angle of grating rotation 
measured from the bisector of the angle 2 K and has the same sign as the spectral order m. 
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and 0, and 0 2 are related to Xj and X 2 through equation (13). This is equivalent to imposing on 
Fj j K the condition 


$ 2 — r 

lijk = / Fjjk d0 = minimum 


( 17 ) 


Determination of Instrumental Parameters p, p', K, q and A 200 

The optimum instrumental parameters are determined in such a way that equation (17) is satisfied 
for the horizontal focusing condition F 200 . This is accomplished by solving the following equa- 
tions simultaneously, 


dI?oo _ - dl 200 _ 

3/3 ~ ’ bp' 

(18) 

9][ioo = q dl 200 _ ^ 

3K " ’ 3A 20 o 

With the values of p", p* and K determined, we proceed to determine the value of q, where q =-—. 

P has been derived by Greiner and Schaffer (ref 4) as 


2 cos K (sin ( d 2 ~0i)) 
(P + P > ) (0 2 -0,) 


(18a) 


where K is the half angle between the entrance and exit slits relative to the vertex of the grating 
0 2 and are the limits of the angular rotation of the grating and 


__ R R 

P andp =~. 

r o to 


q is a very important parameter in that it determines the optimum geometry of the toroid grating 
for the desired wavelength range. Since R is a selected value and the optimum q has been deter- 
mined by (18a), then the vertical radius of curvature of the toroid grating may be calculated using 
equation 


q = 


l 

P 


R_ 

Px 


(18b) 


When the optimum values of p, p\ K and q are determined, then each one of the integrals I jj k 
with i + j + k > 2 becomes a function only of Ajj k . Those values of Ajj k ’s which satisfy equation 
( 1 7) are calculated from equation 


9 


( 19 ) 


where 


dljjk =2 f 2 

^ A ijk 9 1 


F ^ F ijk 
F ijk • r ~ d0 = 0 
dA ijk 


A ijk = [B ijk (0 2 )-B ijlc (0,)]/[2(0, -0 2 )- (sin 20 2 -sin 20, 0)] cosK (20) 

and 

B 20 o (0) = (p/3[cos 0 + cos(0 + K) cos (20 + K) 

+ 3 sin (0 + K) sin K] 

(21) 

+ CP' 13) [cos 0 + cos (0 - K) cos (20 - K) 

- 3 sin (0 - K) sin K] - cos 20 cos K 

Bo 2 o (0) = 2 (p + p') cos 0 - q cos 2 0 cos K (22) 

B 300 = (p/4 ) 2 [sin [40 + 3K] -4 cos 2 (0 + K) sin K - 40 cos K] 

+ (p/4) 2 [sin [40 - 3K] + 4 cos 2 (0 - K) sin K - 40 cos K] (23) 

+ (p/6) [3 sin (0 + 2K] - sin (30 + 2K)J 
+ (p/6) [3 sin (0 - 2K) - sin (30 - 2K)] 

B )20 (0) = (’/2p)2 [sin (20 + K) - 20 cosK] 

+ fl/-.75V rein fOfl - IO - Oft me V] 

(24) 

+ — [3 sin (0 + 2K) - sin (30 + 2K)] 

6 

qp’ 

+ — [3 sin (0 - 2K) - sin (30 - 2K)] 

6 

and so on. 

Determination of the Recording Parameters p c , p d , 7 and 8 and Minimization of Astigmatism 
and Coma 

The recording parents depend upon the aberrations to be minimized. With a predetermined effec- 
tive grating constant and wavelength range, the following equations are solved simultaneously for 
minimizing astigmatism and one coma-type aberration. 
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( 25 ) 


sin 6 — sin 7 = 


Xo 


fioo (Pc.T)-fioo (Pd, 5 ) = (W 0 ) A 20 o 
fo JO (Pc, q, 7 ) - *020 (Pd, q, 5 > * (^ 0 /«t) Aojo 
( 300 (Pc, 7 ) “ (30 (Pd , 5 ) = ( Xo /») A300 
and for minimizing coma-type aberrations 


sin S - sin 7= — 
a 


fj 00 (Pc, 7 ) (200 (P D , 6 ) = (* 0 /«) Ajoo 


f 3 0 0 (Pc , 7 )-f 3 00 (pD»®)~(*°/a) A3 00 

fl 20 (Pc,q,7)~flJ0 (PD,q, 8 ) = ( XO /o) A, 20 


( 26 ) 


Before making any attempt to solve equations ( 25 ) or ( 26 ) it is necessary to investigate the condi 
tions under which equations ( 25 ) or ( 26 ) can have real solutions forp c and p D because H 300 
and H, 20 are quadratic in p c and Pd and their values depend onA 300 and Ai 2 o and therefore 
on the values of 02 in I 30 o and (no- The condition may be stated in such a way that two quad- 
ratic equations of p c or Pd resulting from equations ( 25 ) and ( 26 ) should not have imaginary 
roots. To fulfill this condition, A 300 and A, ao must satisfy the equation 


b 1 — 4ac ^ O 


(27) 


where for A ao o 

a = (sin 8 cos 2 7 - sin 7 cos 2 6) cos 2 7 

b = [2p cos 7 sin 8 + sin (5 - 7) cos 6] cos 7 ( 28 ) 

c = p(p + cos 7) sin 6 + A 3 oo cos 2 6 (sin 8 - sin 7) 
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and for A! 20 


a = cos 4 7 sin 5 - sin y cos 4 8 
b = [ 2 p cos 7 sin 5 + sin (5 - 7 ) cos 3 8 ] cos 7 

(29) 


c = p (p sin 5 + cos 3 8 ) sin 5 + A! 20 cos 4 5 (sin 5 - sin 7 ) 
p = A 200 (sin 8 - sin 7 ) + cos 7 - cos 8 . 

when condition (27) is satisified, then p c and pd are solutions to equations (25) and (26). 

Up to this point, we have been describing the general method of designing aberration corrected 
holographic toroid gratings for Seya-Namioka type monochromators which are not interchange- 
able with conventional toroidal gratings. We will now deal with the design method where both 
types of gratings are interchangeable. 

Modified Method 


The modified method of designing holographic toroidal gratings for the Seya-Namioka mono- 
chromator has a practical advantage in that it may now be designed so that it is interchangeable 
with the conventionally ruled toroidal grating having the same groove frequency and radius of 
curvatures R and p t . This is possible provided the same instrumental constants p, p', q and K 
used for the conventionally ruled gratings are used to design the holographic toroidal grating. The 
instrumental constants for the conventional toroidal grating are determined by solving the 


equation 


la 


00 


= / 2 Mloo 
0 1 


d 0 = minimum 


dlaoo _ ~ dlaoo _ q 

bp ' bp ' 


and 



(30) 


The design equations given by equations (17)-(26) must be modified in part in order to accom- 
modate condition (30). The modification required is to replace equation ( 1 8) with equation (30) 
and equation ( 20 ) with (ijk) = ( 200 ) and the rest of the procedure remains the same. 

The design of toroidal holographic gratings and optical systems for Seya-Namioka monochro- 
mators have been presented. Design equations for minimizing aberrations such as astigmatism and 
coma have been derived. By now, you must have observed the the design equations for concave 
toroid gratings arc very similar to the design equations for concave spherical gratings. In fact, if 
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the parameter q has a value of one; ie p t = R, then the toroidal equations are identical to those 
equations for the sphere. The advantage of the concave toroid grating relative to the concave 
spherical grating is that astigmatism can be minimized for larger angles of incidence and at shorter 
wavelengths where under the same geometrical conditions the concave spherical grating cannot. 

Section 2 


Section 1 described the theory of designing holographic concave toroidal gratings that minimize 
aberrations for a selected wavelength range. In this section, we will apply the theory and present 
a computer program that will calculate all the necessary parameters for using and ruling the holo- 
graphic toroid grating for the Seya-Namioka monochromator. 

Program Description: 

The computer program consists of a main program and four subroutines. Fig (4) is a flow chart 
diagram of the program. 

The purpose of the main program is to direct the flow of calculations depending upon which option 
is selected. If option LMN=2 is selected, then the optimum value of the instrumental parameter 
angle 2K will be calculated. However, if the program user chooses to specify an angle 2K as an 
input value, then option LMN=1 is used. 

Two toroidal holographic grating design methods have been developed: the general method and the 
modified method. If the general method is chosen then the option called I ETA= I is used. When 
using this option, all instrumental parameters are related only to the toroidal holographic grating. 
However, if the modified method is selected, then the option IETA=0 is used in the program. 

When using this option, all instrumental parameters calculated are the same for both the holo- 
graphic and conventional gratings. This means that the holographic and conventional toroid 
gratings are interchangeable if they have the same groove frequency and radii of curvatures. (R 
and p t ). The difference is that the holographic toroid grating minimizes aberrations and the con- 
ventional toroid grating only reduces astigmatism. When the design method has been selected and 
the option has been decided upon, the following steps of calculation are performed: 

Step A: Read input cards 

Step B: Gall INSTRM to calculate the optimum angle 2K, p, p' and p t . Then calculate Bjj k and 
Ajj k from equation (20) where LMN=2 or calculate the functions Bjj k and Ajj k using the selected 
input value for 2K when LMN=1. 

Step C: Call ASTGMS to calculate the recording parameters P c , Pd » 7 and 5 to minimize astig- 
matism ijk = (020) and one coma term tjk = (300). 

Step D: Call COMA to calculate the recording parameters p c , Pd . 7 and 5 to minimize coma 
aberrations ijk = (300) and ijk = (120). 
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Step E: Call PRFORM to display the numerical values of each aberration of the holographic toroid 
grating and compare with that of the equivalent conventional grating as a function of wavelength. 
Note that the program automatically calculates the recording parameters for both the minimiza- 
tion of astigmatism and coma. This feature was included in the program so that the program user 
may select which aberration minimization will best produce the optimum minimization for all 
aberrations through out the desired wavelength range. 

The procedure stated above is supported by four subroutines. 

Subroutine INSTRM 

Subroutine INSTRM determines the optimum values of the instrumental parameters p, j o', 2K, and 
a 200 by simultaneous solution of equation ( 1 8). A modified form of Newton’s method is used to 
obtain the solution. With a knowledge of p, p 1 and 2K the optimum vertical radius of curvature 
p, for the toroid is calculated. After all the optimum instrumental parameters are determined, 
then the values of the B S j k are calculated which in turn are used in equation (20) to determine 
the values of A^ ’s. Note that if a prescribed input value of 2K were used, the selection of option 
LMN=1 would by-pass the optimization procedures when calculating the other instrumental param- 
eters. 

Subroutine ASTGMS 

The purpose of this subroutine is to obtain the recording parameters such that astigmatism (ijk) 

= (020) and one coma term (ijk) = (300) are minimized over the given wavelength range X, < 

X < X 2 . This is achieved by simultaneous solution of equation (25). Before proceeding with the 
solution, it is first verified that the value of NXo is less than the maximum value of |(sin 5 - sin 
7)1, namely 2. If not, an error message is printed stating “error, GAMA out of bounds.” Now, 
if NXo is less than 2 then the values of a, b, and c of equation (28) or (29) are computed and the 
discriminate is checked to ensure that b 2 - 4ac > 0. The solutions of the quadratic equations are 
computed and checked to make sure tha t the positive v alue has been used, because the calculations 
show that one of the two roots [(- b ±y^(b 2 - 4ac)/2al is always positive and the other is always 
negative. When the discriminate has been verified to be positive, then the optimum values of 
Pc Pd > 7 and 8 are determined by the iterated solution of equation (25) using Newton’s method. 

Subroutine COMA 

The purpose of this subroutine is to obtain the recording parameters such that the coma terms 
ijk = (300) and yk = ( 1 20) are minimized for a given wavelength range X, < X < X 2 . This is 
achieved by a simultaneous solution of equation (26). We use the same procedure used in sub- 
routine ASTGMS to verify that the discriminant is positive and then proceed to determine the 
optimum values p c , Pd » 7, and 6 by the iterated solution of equation (26) using Newton’s 
method. 
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MAIN PROGRAM FLOW DIAGRAM 



INPUT 

i < AX, C, N, R, K, 
I ETA, LMN, QTI, QTF, 
DOT, IQ 
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Subroutine PRFROM 


This subroutine calculates the performance of the holographic toroid grating using the calculated 
instrumental and recording parameters and then displays the values of each aberration as a function 
of wavelength in the form of a table. The performance of a conventional (mechanically ruled) 
grating with the same groove frequency is also shown in the table for comparison. 

Section 3. DATA SET AND DESCRIPTION OF DATA CARDS 

This section discusses the input data required for a successful operation of the program. 

(A) Data Set 

The following parameters are required as input for the program. 

1. X Q = wavelength of recording laser light in angstroms. 

2. X L = lower wavelength limit in angstroms. 

3. X u = upper wavelength limit in angstroms. 

4. AX = wavelength interval in angstroms which aberrations will be displayed. 

5. C = angle 2K (0.0 is used with option LMN=2). 

6. N = groove frequency (number of lines per mm). 

7. R = Radius of curvature in mm. 

8. m = order of diffraction where m = ± 1 , ± 2, etc. 

9. 1ETA: Parameter specifying the grating design method. IETA = 1 specifies general method. 

IETA=0 specifies the modified method. 

1 0. QT1 = Initial value of the torus parameter qj. 

1 2. QFT = Final value of torus parameter q f . 

1 3. DQT = Interval at which the torus parameter Aq will change. 

14. IQ = Parameter to specify the torus radius options. 

IQ=2; Torus parameter q is selected by user. 

IQ= 1 ; Toroid radius p, in mm is automatically computed 
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Note: If IQ=1 is selected, then QTI, QTF and DQT are NOT punched on the data card. The 

computer automatically bypasses these parameters and calculates the optimum q values and 
p t . However the value 1 for IQ must be punched in the appropriate column on the card. 

If IQ=2 is selected, then QTI, QTF, and DQT must be punched on the data card in the 
correct order. For this case, the computer will calculate the p, based on the values selected 
for the range of q’s desired. 

(B) Description of Data Cards and Sample Data Cards 

The first seven input parameters must be punched on the data cards in the order shown in (A) 
starting in column 1 in FI 0.1 format with ten columns reserved for each of the seven parameters 
and each parameter must include a decimal. The next three parameters are in 12 format starting in 
column 71. The parameters QTI, QTF, DQT, and IQ are punched on the second card in format 
F10.0 and 12 respectively. Fig (5) shows a typical two card data set to illustrate the above 
description. 

Note: ( 1) The program user may include as many data cards for a single run as desired. 

(2) The values of all 14 parameters must be punched for each additional set. 

Section 4 


This section presents a typical example for using the program and also describes the output data in 
detail. 

Example: Determine the optimum instrumental and recording parameters for designing a concave 
toroid holographic grating for a Seya-Namioka monochromator that minimizes astig- 
matism or coma and has the following specifications: 

1. Interchangeable with conventional concave toroid grating. 

2. Groove frequency = 550 L/mm 

3. R = 1,000 mm 

4. Operating wavelength range = 0.0A to 1600A. 

5. Order = -l 

6. C= 142° 

Solution: 


( 1) indicates that the modified method is used; IETA=0. (6) says that option 
LMN=1 is used. 
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Figure 5. Typical Data Card Set 
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Jhe first data card is punched using data supplied above and punched in the order listed below: 
X 0 = 4579.3 
X L = 0.0 
X„ = 1600.0 
AX = 100.0 
C= 142.0 
N = 550.0 
R= 1000.0 
K = -1 
IETA = 0 
LMN = 1 

The following is punched on the second data card: 

QIT = 1.0 
DQT = 1.0 
QTF = 2.0 
IQ= 1 

The output data for the example cited is shown in fig 6 and described as follows: 

Line 1 Shows the input parameters; order = -1 is the order of the grating, N = groove frequency 
in grooves per mm, R = horizontal radius of curvature of the grating in mm. IETA = 
identifies the design method and LMN = 1 shows that the angle C has been specified. 

Line 2 shows the input values for X 0 = LAMDA 0 (Laser Wavelength), X L = LAMDA 1 
(lower wavelength limit), X u = LAMDA 2 (upper wavelength limit). 

Line 3 is self-explanatory. 

Line 4 shows the calculated instrumental parameters; C = 2K, p = RHOA, p' = RHOB, 7 0 = SMLR 
in mm, and r^, = SMLR in mm. These parameters are defined in the text. 
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Line 5 shows the values of the torus parameters: Q = q value defined in the text, p t = value of 
- the vertical toroid radius. 

Line 7 is self-explanatory and refers to line 8. 

Line 8 shows the optimum recording parameters; y = GAMA, 5 = DLTA, p c = RHOC, and p D = 
RHOD. These parameters are defined in the text. 

Line 9 shows the titles for the table showing the performance of the grating as a function of wave- 
length. The table also shows the comparison of the holographic toroid grating with the 
toroid conventional grating. 

\ = LAMDA is the wavelength at which the grating is evaluated. 

F 20 o X 10 -4 is the horizontal focus. 

F 0 2 0 is the vertical focus and astigmatism. 

F 300 is the first coma term. 

F, 20 is the second coma term. 

F 400 is the higher order focusing term. 

F 220 >s the toric aberration. 

F 040 is the toric aberration of a higher order 

Line 1 0 shows the holographic grating performance (Modified Method), the value of the wave- 
length and ine values of the aberrations. 

Line 1 1 shows the performance of the conventionally ruled grating (MECH. Ruled). 

Line 12 shows that the following information refers to minimizing coma type aberrations for the 
same instrumental parameters shown in line 4. 

Line 13 is self explanatory and refers to line 14. 

Line 14 shows the recording parameters y - GAMA, 5 = DLTA, p c = RHOC, and p D = RHOD 
required 1 to fabricate a toroid grating to minimize coma. 

Lines 15, 16, and 17 are defined in the same way as lines 9, 10, and 1 1 above. The values in this 
table refer to a grating designed to minimize coma. 
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Figure 6. 
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